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Abstract 

Using the Batalin-Vilkovisky formalism, we study the Ward identities and the equations of 
gauge dependence in potentially anomalous general gauge theories, renormalizable or not. A cru¬ 
cial new term, absent in manifestly nonanomalous theories, is responsible for interesting effects. 
We prove that gauge invariance always implies gauge independence, which in turn ensures per¬ 
turbative unitarity. Precisely, we consider potentially anomalous theories that are actually free of 
gauge anomalies thanks to the Adler-Bardeen theorem. We show that when we make a canonical 
transformation on the tree-level action, it is always possible to re-renormalize the divergences 
and re-fine-tune the finite local counterterms, so that the renormalized P functional of the trans¬ 
formed theory is also free of gauge anomalies, and is related to the renormalized P functional of 
the starting theory by a canonical transformation. An unexpected consequence of our results is 
that the beta functions of the couplings may depend on the gauge-fixing parameters, although 
the physical quantities remain gauge independent. We discuss nontrivial checks of high-order 
calculations based on gauge independence and determine how powerful they are. 
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1 Introduction 


The Ward-Takahashi [1, 2] and Slavnov-Taylor [3, 4] identities are relations among the correlation 
functions of quantum field theory, and follow from gauge and global symmetries. They are usually 
studied in theories that are manifestly nonanomalous, that is to say admit a manifestly gauge 
invariant regularization technique, for example QED and nonchiral Yang-Mills theories. Chiral 
gauge theories, such as the standard model, are potentially anomalous, because they do not admit 
a manifestly gauge invariant regularization technique. The Adler-Bardeen (AB) theorem [5, 6, 7] 
is the main tool that can establish whether a potentially anomalous theory is in the end truly 
anomalous or nonanomalous. It ensures that, if the gauge anomalies are trivial at one loop, they 
can be cancelled to all orders. 

The potentially anomalous theories that are actually free of gauge anomalies thanks to the 
Adler-Bardeen theorem will be called AB nonanomalous. In this paper, we study the Ward 
identities of the AB nonanomalous general gauge theories, including the nonrenormalizable ones, 
and clarify the relation between gauge invariance and gauge independence. Our investigation 
upgrades the ones available in the literature in several respects. 

Gauge invariance and gauge independence are two different concepts, to the extent that a 
functional can be gauge invariant and gauge dependent at the same time. For example, the renor¬ 
malized action of non-Abelian Yang-Mills theory contains a term propotional to Zj^ f 
where is the field strength and Za is the wave function renormalization constant of the gauge 
field. This expression is gauge invariant, but not gauge independent, because Za may depend on 
the gauge-fixing parameters. 

Yet, the two concepts are related to each other, and crucial to prove perturbative unitarity. 
Gauge invariance is necessary, because its violation makes unphysical degrees of freedom, such 
as the longitudinal photons, propagate. On the other hand, gauge independence is important, 
because it allows us to switch back and forth between gauges that exhibit perturbative unitarity, 
but do not have good power-counting behaviors (such as the Goulomb gauge), and gauges that 
have good power-counting behaviors, but do not exhibit unitarity (such as the Lorenz gauge). 
The Lorenz gauges are very convenient to make calculations and prove theorems to all orders. 
They make renormalizability manifest, when the theory is power-counting renormalizable. When 
the theory is nonrenormalizable, they make the locality of counterterms manifest. However, the 
Lorenz gauges hide unitarity, because they introduce unphysical, propagating degrees of freedom, 
such as the longitudinal components of the gauge fields and the Fadeev-Popov ghosts. This is 
where gauge independence plays a key role, because it ensures that every physical quantity can be 
equivalently defined by using the Goulomb gauge, where the propagators have no unphysical poles 
and perturbative unitarity is manifest. The equivalence of the two gauges allows us to loosely say 
that “the unphysical degrees of freedom of the Lorenz gauges compensate one another and drop 


2 



out of the physical quantities”. 

Thus, in quantum field theory we need both gauge invariance and gauge independence. If a 
theory is AB nonanomalous, it is by definition gauge invariant. It is not obvious that the Adler- 
Bardeen theorem also ensures that the physical quantities are ultimately gauge independent. 
Is it so, or do we need extra assumptions to ensure that the physics does not depend on the 
gauge fixing? Among other things, in this paper we answer this question by proving that gauge 
invariance always implies gauge independence. 

In our approach, the Ward identities of AB nonanomalous general gauge theories are corrected 
by a term that is absent in manifestly nonanomalous theories. The correction is evanescent at 
the bare level, but can generate finite corrections at the renormalized level, by simplifying some 
divergences. One of the main consequences is that the beta functions of the couplings can depend 
on the parameters introduced by means of the gauge fixing. However, the physical quantities are 
protected, and remain gauge independent. 

We study how the renormalized T functional T^ depends on the parameters introduced by 
the canonical transformations of fields and sources. Canonical transformations encode field re¬ 
definitions and changes of the gauge fixing, both of which are expected to have no effect on the 
physical quantities. When we speak of “gauge dependence” we refer to the dependence on all 
types of parameters introduced by a canonical transformation, including those associated with 
field redefinitions. 

We work out how a canonical transformation on the (bare) action S affects the renormalized 
r functional T^^. After the transformation, the theory must be renormalized anew. We show that 
in this process of re-renormalization, it is always possible to redefine the subtraction scheme, by 
fine-tuning the finite local counterterms, so that the transformed theory is also AB nonanomalous. 
Moreover, the gauge dependence of the transformed T/j is encoded into a canonical transformation, 
up to evanescent corrections. 

This result allows us to prove that the physical quantities are gauge independent. However, 
quantities that are useful for intermediate purposes, such as the beta functions of the couplings, 
are normally gauge dependent. Their gauge dependence can be absorbed inside finite redefinitions 
of the couplings. 

In manifestly nonanomalous theories we are, to a large extent, free to use a preferred subtrac¬ 
tion scheme, such as the minimal one, both before and after the canonical transformation. The 
physical quantities and the beta functions of the couplings are unaffected by the transformation 
(see for example [8]). In AB nonanomalous theories, instead, we can use a preferred subtraction 
scheme neither before, nor after the transformation. Before the transformation, we need to choose 
a specific class of subtraction schemes to take advantage of the Adler-Bardeen theorem and cancel 
the gauge anomalies to all orders. After the transformation, we need to choose (another) specific 
class of subtraction schemes, to enforce the cancellation of gauge anomalies again. In this process. 
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some gauge-fixing parameters move out of the gauge-fixing sector into another unphysical sector, 
the one encoded by the choice of the subtraction scheme. The result is that the beta functions 
are gauge dependent, in general. Nevertheless, we can make their gauge dependences disappear, 
if we specify the new subtraction scheme even further. 

Both gauge invariance and gauge independence can be used to make powerful checks of high- 
order calculations. As said, a consequence of our investigation is that in AB nonanomalous 
theories, including the standard model, the beta functions of the couplings are not completely 
gauge independent. We show that, in spite of this, sufficiently powerful checks of high-order 
calculations are still available. The reason is that the gauge dependence cannot be arbitrary, 
because it cannot affect the physical quantities. 

To keep track of gauge invariance through renormalization, we use the Batalin-Vilkovisky 
(BV) formalism [9]. The gauge invariant regularization techniques commonly used for manifestly 
nonanomalous theories are also convenient to treat AB nonanomalous theories, because they min¬ 
imize the number of terms that are potentially anomalous. In this paper we use the dimensional 
regularization [10], or any regularization technique that underlies the dimensional one, such as the 
chiral dimensional (CD) regularization of ref. [11] and the (chiral)dimensional/higher-derivative 
regularization of refs. [6, 7, 11, 12], obtained by merging the (chiral) dimensional one with the 
covariant higher-derivative regularization of ref. [13]. We recall that the CD regularization is 
particularly convenient for studying nonrenormalizable theories, to avoid certain ambiguities that 
show up when we extract the divergent parts of the BV antiparentheses (X, Y) of two functionals 
X and y, as well as other nuisances that the ordinary dimensional regularization is responsible 
for. 

We also take the chance to revisit some known issues under our perspective. 

Before presenting our results in more detail, we comment on the existing literature on related 
subjects, and explain the upgrades we make. Most studies of gauge dependence have been focused 
on renormalizable theories [14], or nonrenormalizable, but nonchiral, theories [15, 16], where the 
problem is much simpler (see appendix E). We want to develop an approach that also applies 
to nonrenormalizable chiral theories, to include the standard model coupled to quantum gravity. 
In our opinion, it is not necessary to wait for the ultimate theory of quantum gravity to prove 
general statements about it. The other investigations of gauge dependence we are aware of use the 
so-called algebraic approach to renormalization [17]. The main feature of the algebraic approach 
is that it does not make use of an explicit regularization technique. Instead, it relies on tools such 
as the “quantum action principle” [18]. 

We think that it is important to develop more standard approaches to the problem of gauge de¬ 
pendence, like the one of the present paper, which uses the dimensional regularization or modified 
versions of it. For example, anomalies have taught us that working without an explicit regulariza¬ 
tion may not be completely safe. Another advantage of using an explicit regularization is that we 
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can identify convenient subtraction schemes, where simplihcations occur and several properties 
are easier to deal with, to all orders in the perturbative expansion. Examples are those provided 
by refs. [6, 7], where it was shown that in suitable subtraction schemes the gauge anomalies 
automatically vanish from two loops onwards, if they cancel out at one loop. By construction, it 
is not possible to identify special subtraction schemes in regularization-independent approaches. 

Now we state the main results of our investigation. We study canonical transformations that 
are continuously connected with the identity. Their generating functionals have the form 

F{^,K',9) = J + (1.1) 

where 6 denotes the “gauge parameters”, which are associated with both changes of field variables 
and changes of the gauge fixing. In the first part of our analysis, we prove the main theorem, 
which states that if the theory is AB nonanomalous at 0 = 0, after making the canonical trans¬ 
formation (1.1) it is always possible to re-renormalize the divergences and re-fine-tune the finite 
local counterterms, continuously in 9, so that the equations 

(ABT) (rKe,r^;e) = 0(e), (1.2) 

(GDE) ^_(r^,,(Q^,)) = 0(e) (1.3) 

hold for arbitrary 9, where Tjiq is the renormalized E functional of the transformed theory and Qjiq 
is a suitable renormalized local functional. The right-hand sides of both equations are (generically 
nonlocal) functionals that vanish when the continued spacetime dimension D = d — e tends to the 
physical spacetime dimension d. We denote such functionals by 0(e) and call them “evanescent”. 

Equation (1.2) ensures that the theory is AB nonanomalous for arbitrary values of the gauge 
parameters. Thus, it encodes gauge invariance. Eormula (1.3) is the equation of gauge dependence, 
and follows from the generalized Ward identities. The equations (GDE) can be integrated to show 
that the entire gauge dependence of T^je can be absorbed inside a (convergent, but generically 
nonlocal) canonical transformation, up to 0(e). The results encoded in formulas (ABT) and 
(GDE) are so general that they do not require any particular assumption (see section 3). 

We also derive the equations of gauge dependence at the level of the renormalized action and 
show that RG invariance is preserved by the canonical transformation. 

A simple, but important application of the theorem is to power-counting renormalizable chiral 
gauge theories gauge-fixed by means of a nonrenormalizable gauge fixing. We show that the theory 
remains renormalizable in a nonmanifest form, because the parameters of negative dimensions 
introduced by the gauge fixing do not propagate into the physical sector. Another application of 
the theorem is a crucial step in the proof of the Adler-Bardeen theorem for nonrenormalizable 
theories [7]. 
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In some situations, we can prove formula (ABT) for arbitrary values of a certain gauge pa¬ 
rameter 0 within a given class of subtraction schemes. Then, it is not necessary to re-renormalize 
the divergences and the re-fine-tune the finite local counterterms. Under the assumption that the 
theory satisfies a certain cohomological property, which is a generalized version of the well-known 
Kluberg-Stern-Zuber conjecture [19], we can derive an more specific version of equations (GDE), 
which reads 

(GDE2) ^ - (T^,, = 0(e), (1.4) 

i 

where A* are the independent parameters of the classical action, pi are constants that depend on 
Aj and the other parameters of the theory and Hp^Q is a renormalized local functional. 

We can write (1.4) in the form (1.3) by suitably “evolving the parameters A in the 9 direction”. 
Such redefinitions encode how the beta functions of the couplings depend on 9. 

So far, the Adler-Bardeen theorem has been proved in a variety of cases. The original 
proof given by Adler and Bardeen [5] was designed to work in QED. Most generalizations to 
renormalizable non-Abelian gauge theories used arguments based on the renormalization group 
[20, 21, 22, 23], which work well unless the hrst coefficients of the beta functions satisfy pecu¬ 
liar conditions [23] (for example, they should not vanish). Then there exist algebraic/geometric 
derivations [24] based on the Wess-Zumino consistency conditions [25] and the quantization of 
the Wess-Zumino-Witten action. Another method to prove the Adler-Bardeen theorem in renor- 
malizable theories is obtained by extending the coupling constants to spacetime-dependent fields 
[26]. A proof that covers all power-counting renormalizable gauge theories was given in ref. [6]. It 
was obtained by elaborating on a previous proof [12] given for quantum field theories that violate 
Lorentz symmetry at high energies (in particular, Lorentz violating extensions of the standard 
model) and are renormalizable by weighted power counting [27]. Recently, the proof of [6] was 
further extended in ref. [7], to include a large class of nonrenormalizable theories, such as the 
standard model coupled to quantum gravity. We emphasize that a byproduct of our investigation 
is that the standard model, coupled to quantum gravity or not, is perturbatively unitary, and so 
are most of its extensions. 

The paper is organized as follows. In section 2 we compare the Ward identities of chiral 
and nonchiral gauge theories, and illustrate the crucial new term that appears when the theory 
is potentially anomalous. In section 3 we prove the main theorem of this paper, by deriving 
and integrating the equations of gauge dependence in AB nonanomalous theories. We show that 
every canonical transformation on the classical action is mapped into a canonical transformation 
on the renormalized T functional, provided that the finite local counterterms are appropriately 
re-fine-tuned. We also integrate the equations of gauge dependence. In section 4 we derive the 
equations of gauge dependence of the renormalized action. In section 5 we prove that the canonical 
transformation preserves RG invariance and discuss two applications of the main theorem. In 
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section 6 we study the gauge dependence of the beta functions in detail. In section 7 we explain 
how to switch off the ghosts, the antighosts, the Lagrange multipliers for the gauge fixing, and 
the sources for the symmetry transformations, and get to the physical quantities, collected into 
a “physical” L functional Lph- We derive the (nonlocal) gauge symmetry of Lph, and prove that 
it closes off shell. Finally, we prove that Fph is gauge independent, up to field redefinitions, and 
perturbatively unitary. In section 8 we investigate the checks of high-order calculations provided 
by gauge independence and estimate how powerful they are. Section 9 contains our conclusions. 
In the appendices we prove some properties used in the paper, recall earlier results and collect 
some reference formulas for the standard model coupled to quantum gravity. Moreover, we revisit 
the gauge dependence of manifestly nonanomalous theories in the light of the new results. 


2 Generalized Ward identities 


In this section we fix some notation, recall the main properties of the Batalin-Vilkovisky formalism 
for general gauge theories [9] and derive the generalized Ward identities. 

Let D = d — £ denote the continued, complex dimension of spacetime, and d the physical 
spacetime dimension. The ZD-dimensional spacetime manifold is split into the product 
of the ordinary d-dimensional spacetime R*^ times a residual (—e)-dimensional evanescent space, 
R“^. The spacetime indices ... of vectors and tensors are split into the bar indices fi,!?,.. 
which take the values of 0 , 1 , • • • ,d— 1 , and the formal hat indices fi,!), ■ ■., which denote the 
R“^ components. For example, the momenta are split into the pairs pf^, also written as 
p^, p^, and the coordinates are split into x^, x^. The formal flat-space metric is split 
into the usual d x d flat-space metric =diag(l, — 1 , • • • , — 1) and the formal evanescent metric 
Pup = —5pp. The off-diagonal components ppp vanish. The evanescent components are contracted 
among themselves by means of the metric ppp, so for example p^ = p^pppp'^■ Full SO{l,D — 1) 
invariance is lost in most expressions, replaced by SO{l,d — 1 ) x SO{—£) invariance. 

We recall that in the CD regularization the fields have strictly d-dimensional components. 
The metric tensor is block-diagonal: the diagonal blocks are gpp{x) and ppp, while gpp = 0. 
Moreover, the 7 matrices are strictly d dimensional, and satisfy the usual Dirac algebra { 7 “, 7 ^} = 
27 “^, where the indices d,b,... refer to the Lorentz group. If d = 2A: is even, the d-dimensional 
generalization of 75 is defined as 


7 = 


and satisfies 7 ^ = 7 ) 7 ^ = 1 - The left and right projectors Pl = (1 — 7 )/ 2 , Pr = (1 -|- j )/2 are 
defined as usual. The tensor and the charge-conjugation matrix C also coincide with the 

usual ones. 
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The set of fields = {(p’’, C, C, B} contains the classical helds cp, the Fadeev-Popov ghosts 
C, the antighosts C and the Lagrange multipliers B for the gauge fixing. An external source 
with opposite statistics is associated with each and coupled to the transformations i?"(d>). 
If X and Y are functionals of and K, their antiparentheses are defined as 


Y)- f 


( 2 . 1 ) 


where the integral is over spacetime points associated with repeated indices and the subscripts I 
and r in 5i and 6r denote the left and right functional derivatives, respectively. 

The action S should solve the master equation (S, S) = 0 in D dimensions, with the “boundary 
condition” S{^,K) = Sc{(p) SitC = C = B = K = 0, where Sc{<p) is the classical action. 

If the gauge algebra closes off shell, there exists a choice of field/source variables such that 
the non-gauge-hxed solution Sd{^,K) of the master equation has the form 


5d(d>, K) = SM + Sk, Sk{^, K) = -J 


( 2 . 2 ) 


In this case, (Sd, Sd) = 0 splits into the two identities 


/ 




SiScjrP) 


= 0 , 





= 0 , 


which express the gauge invariance of the classical action and the closure of the algebra, respec¬ 
tively. The gauge-hxed solution Sd{^,K) of the master equation reads 


Sdi^, K) = Sc{cP) + {Sk, 'i’) + SK = Sd + {Sk, 'L), (2.3) 


where '!'(<!>) is the gauge fermion, that is to say a local functional of ghost number —1 that encodes 
the gauge hxing. Reference formulas for Sc, Sk and T in the case of the standard model coupled 
to quantum gravity can be found in appendix D. Typically, 'h has the form 

T(<h) = IC (^G{cP, 0 + ^P{cP, e, d)B^ , (2.4) 

where G{(p, f) is the gauge-hxing function, P is an operator that may contain derivatives acting 
on B, and are gauge-hxing parameters. For example, G{(p) = for the Lorenz gauge in 

Yang-Mills theories. Clearly, Sd also solves the master equation {Sd, Sd) = 0 in D dimensions. 

If the gauge algebra does not close off shell, Sd{^,K) is not linear in K and Sd is obtained 
from Sd by applying the canonical transformation generated by 

F{<^,K') = J <^<^K'c + ^{<^>). (2.5) 

In manifestly nonanomalous theories we can solve {S, S) = 0 in Z) dimensions at the regularized 
level. Typically, the solution coincides with (2.3). In potentially anomalous theories, instead, we 




cannot achieve this goal. There, the functional Sd{^, K) does solve {Sd-, Sd) = 0 in D dimensions, 
but is not well regularized. The most common reason is the presence of chiral fermions. We can 
deform Sd into a well-regularized action 


S{i^,K) = Sd + S,^ 


( 2 . 6 ) 


by adding an evanescent part S^v that collects suitable regularizing terms [11]. The deformed 
action S does not solve (S', S') = 0 in H dimensions. Instead, it solves the deformed master 
equation 

(5,5) = 0(e), (2.7) 


where the right-hand side denotes terms that vanish for D ^ d. 

Given a generic action 5(<h,i0), the generating functionals Z and W of the (connected) cor¬ 
relation functions are defined by the formulas 


Z{J,K)= + i j =expiW{J,K), 


( 2 . 8 ) 


and the generating functional r(<h,il') = W{J,K) — f of the one-particle irreducible dia¬ 

grams is the Legendre transform of W{J, K) with respect to J. The anomaly functional is defined 
as 


yi = (r,r) = ((5,5)) 


(2.9) 


and collects the set of one-particle irreducible correlation functions that contain one insertion of 

(5,5), where {■■■) denotes the average defined by 5 at arbitrary J . The last equality of (2.9) 
can be proved by making the change of variables 

$“^^>“ + tz7(5,$“) = $“-ro-^, (2.10) 

in the functional integral (2.8), where ro is a constant anticommuting parameter. For the detailed 
proof, see for example the appendices of refs. [6, 8]. See also appendix A. 

Let us explain the meaning of formula (2.9). The functional (5,5) represents the symmetry 
violation, so it is basically the integral of the divergence of the gauge current multiplied by 
the ghosts: 

(5,5) ~2y dPxC{x)^f,J^^{x), 

where the sign means that the right-hand side is written up to terms proportional to the field 
equations and other terms that we can neglect in the present discussion. As said, formula (2.9) 
collects the one-particle irreducible diagrams that contain one insertion of (5, 5) and arbitrary 
external and K legs. The key diagram of this type in four dimensions is the one-loop triangle 
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diagram that is responsible for the well-known ABJ anomaly [28], which arises by considering one 
{S, S) insertion and two external gauge field legs. Amputating those legs, we get 

^{iS,S).Mx)My)) ^ I d^zC{z){d,JP{z)J^{x)My)). ( 2 . 11 ) 

The sign comes from the leg amputation and the fact that we have taken the ghosts out of 
the average, because this is the only way to get nontrivial contributions to anomalies at one loop. 
See ref. [11] for the calculation of the one-loop triangle anomaly in chiral Yang-Mills theories with 
formula (2.9) and the CD regularization technique. 

The Adler-Bardeen theorem is the statement that if the gauge anomalies are trivial at one 
loop, there exists a class of subtraction schemes where they vanish to all orders, that is to say 

Ar = {Tr, Fr) = {{Sr, Sr)) = 0(e), (2.12) 

Sr and T/j being the renormalized action and the renormalized T functional, respectively. The 
right-hand side of (2.12) vanishes for D —)■ d, which ensures that the renormalized T functional 
is gauge invariant in the physical limit. The AB nonanomalous theories are those that admit 
subtraction schemes where (2.12) holds. 

While the AB identity (2.12) ensures gauge invariance, it does not say much about gauge in¬ 
dependence, which is a different statement, namely the property that a certain class of correlation 
functions (that we call “physical”) do not depend on the gauge fixing. 

One way to study the gauge independence is through Ward identities. We begin by recall¬ 
ing how those identities work in manifestly nonanomalous theories, where the master equation 
{S, S') = 0 is satisfied exactly at the regularized level. Let T(<I>) denote a AT-independent, but oth¬ 
erwise completely arbitrary, product of elementary and local composite fields at distinct points. 
By making the change of field variables (2.10) in the functional integral 

J [dd>]Te*^, 

we find 

y [d$] (A, T) e*^ = 0. (2.13) 

We omit details of the derivation, because the proof of this formula is a particular case of the more 
general proof given below. We just stress that it is crucial to use the master equation (S, S) = 0, 
which implies that S is invariant under the field redefinition (2.10). 

Equation (2.13) is the usual Ward identity. For example, if we take T = C{x) d^A^{y) and 
T = C{x)'ip{y)'ip{z) in QED, we can derive the well-known formula Zf,Z^ = 1 that relates the 
renormalization constants Zg and Z^ of the electric charge and the gauge field [29]. 

In this paper, the average (• • •) denotes the sum of connected diagrams. For example, if X and 
Y are local functionals, we have {XY) = {XY)^^— {Z^){Z), where {XY)^^ includes disconnected 
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diagrams. The subscript 0 in (• • • )o means that the correlation functions are evaluated at J = 0. 
An equivalent form of the identity (2.13) is 


((-5,T))o = 0. 


(2.14) 


If we repeat the argument leading to (2.13) without assuming (S, S) = 0, we get the generalized 
Ward identity that we consider in this paper, which reads 

((5,T))o + l((5,5)T)o = 0. (2.15) 

The extra term on the left-hand side of this formula is going to appear in many other contexts 

and is responsible for the new effects anticipated in the introduction. 

To prove (2.15), express T as the product WXi of AT-independent elementary and local com- 

i 

posite fields Aj. Then, consider the functional integral 

y'[dd>]e*^+^^^*"S (2.16) 


where cjj are arbitrary constants. Under the field redefinition (2.10), the action S and the func¬ 
tionals Xi transform as follows: 


S ^ S-w 


/ 


5rS 6lS 

6Ka 


5 + |(5,5), 


Xi —y Xi — tZ7 


< 5.5 6iXi 

6Ka 


Xi + w{S,Xi). 


In the last step we have used the assumption that Xi depends only on the fields 4>. When we 
make the change of variables (2.10) inside (2.16) and divide by (2.16), we get 


/[d^] {'Lj^iS,Xj)aj + §ro(5,5)) 

The left-hand side of this formula is a sum of connected diagrams. Differentiating it once to the 
right with respect to each iTi, ..., fin and setting (Tj = 0 at the end, we project onto the diagrams 
that have one external Uj leg for each i. So doing, we get precisely formula (2.15). 

When the local functionals Xi of the product T = WXi depend on both and K, and the 

i 

sources J are not set to zero, the generalized Ward identities can be worked out from formula 
(2.9), by deforming the action 5 into 5 -|- Ajcii, where ai are constants, and taking the first 
order in all cjjS. 

In particular, if T is equal to a local functional A, it is easy to show that when the action 5 
is deformed into 5 -|- Act, where a is a constant, the T functional deforms into T -|- (A)cj -|- O(cj^), 

while the average (A) of a local functional A deforms into (A)-|-f(AA)rcj-|-0(cj^), where (])([ Aj)r 

i 

denotes the set of one-particle irreducible diagrams that contain one Ai insertion for each i, Ai 
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being local functionals (details are given in appendix A). Expanding (r,r) = {{S,S)) in powers 
of a and taking the first order of the expansion, we obtain the identity [8] 

((5,A)) + l((5,5)A)r = (r,(A)). (2.17) 

Both sides of (2.17) are viewed as functionals of and K (rather than functionals of J and K). 
Note that, in particular, {X) = (X)r. 

Repeating the derivation for T = XY^ where X and Y are both local functionals, we get the 
identity 

((5, Xy))r + \{{S, S)XY)r = (E, (Ay)r) - i(-l)^- ((X), (y)) + i(-l)^- ((X, X)), (2.18) 

where £x denotes the statistics of the functional X (which is 0 if X is bosonic, 1 if it is fermionic). 
When T is the product of more local functionals, we can proceed similarly. 

An important application of the generalized Ward identities is the derivation of the equations 
of gauge dependence, which tell us how the generating functional T depends on the gauge param¬ 
eters. We first recall such equations in manifestly nonanomalous theories and then switch to AB 
nonanomalous theories. 

In manifestly nonanomalous theories {S, S) = 0 in D dimensions and 5ev = 0, 5 = S^- The 
functional T satisfies the equation 

(dr / (d S' \ 

^ = = = (2.19) 

where f is any gauge-fixing parameter and = d^/df is the ^-derivative of the gauge fermion 
'k. The first equality is obvious. The second equality follows from formula (2.3). Indeed, recalling 
that the parameters f are contained only in 'll, we have dS/df = {Sk, 'k^) = {S, dt^). The third 
equality follows from formula (2.17). 

More generally, if 9 denotes any gauge parameter, introduced by a canonical transformation 
generated by (1.1), we find 

(dr / (d S \ 

^ = (^) = ((^,0.)) = (r,(g.)), (2.20) 

where Qe is the derivative T’(<k, K', 9) with respect to 9, reexpressed as a functional of <k and K. 

Equations (2.19) can be renormalized and integrated (see [8] and appendix C). The result 
is that the ^ dependence can be absorbed into a canonical transformation on F. Therefore, the 
contributions due to the right-hand side of (2.19), which are in general nonvanishing, do not affect 
the physical quantities, for example the S-matrix elements. See subsection 7.3 for details. 

In AB nonanomalous theories the equations of gauge dependence are corrected by an extra 
term, which corresponds to the extra term of (2.15). Formula (2.20) turns into [8] 

% = {%) = ( 2 . 21 ) 


12 



Assuming that the primes denote the ^-independent quantities, the second equality of (2.21) 
follows from formula (A.6) recalled in appendix A, since dS'/dO = 0. The last equality of (2.21) 
follows from formula (2.17). 

The identities (2.15), (2.17), (2.18) and (2.21) are so general that they also hold in truly 
anomalous theories. However, their most interesting applications are to AB nonanomalous theo¬ 
ries, which are the main focus of this paper. 

In the next sections we are going to renormalize the equations (2.21) and integrate their 
renormalized versions. The nontrivial part of this task is to work out the effects of the last 
term of formula (2.21). The result is that the 0 dependence can be absorbed into a canonical 
transformation on the renormalized T functional T^, provided that the finite local counterterms 
are appropriately fine-tuned. 

We stress again that gauge invariance, which is expressed by formula (2.12), does not imply 
gauge independence in an obvious way. However, in this paper we prove that ultimately it does. 
Gauge independence allows us to prove the perturbative unitarity of the theory (see subsection 
7.4). 

Before concluding this section, we make some remarks to emphasize the role played by the 
evanescent terms 0(e) in our discussion. With respect to the limit D ^ d we can distinguish 
divergent, nonevanescent and evanescent terms. A contribution is called “nonevanescent” if it 
has a regular limit for D ^ d and coincides with the value of that limit. In the (ordinary, 
as well as chiral) dimensional regularization the evanescences can be of two types: formal or 
analytic. Analytically evanescent terms are those that factorize at least one e, such as eFfipF^^, 
^tc., where ipL is a left-handed fermion. Formally evanescent terms are those that 
formally disappear when D ^ d, although they do not factorize powers of s, such as The 

divergences are poles in e, and can multiply either nonevanescent terms or formally evanescent 
terms. In the latter case they are called divergent evanescences. An example is It is 

convenient to subtract away the divergent evanescences like any other divergences. 

In most derivations it is necessary to extract the divergent parts of functionals and antiparen¬ 
theses of functionals. We have to take some precautions to ensure that this operation can safely 
cross the antiparentheses, so that for example (S, A)ciiv = (<5, Wdiv)- The first thing to do is 
define the classical action (2.6) so that it does not contain analytically evanescent terms, but only 
nonevanescent and formally evanescent terms, multiplied by e-independent coefficients. In this 
way, S does not contain dangerous e factors that could simplify the divergences of X inside (S', A). 
For the same reason, it is convenient to use the chiral dimensional regularization of [11], instead of 
the ordinary dimensional regularization. In particular, we must use the CD regularization when 
the theory in not power-counting renormalizable. So doing, we avoid a number of ambiguities 
that would complicate our operations. For details on this subject, see refs. [6, 11]. 
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3 The theorem of gauge dependence 

Consider a general gauge theory with action S{^, K,uj), where tc denotes its parameters. Let Sr 
denote the renormalized action and F/j the renormalized F functional. Assume that the theory is 
AB nonanomalous, i.e. 

{TR,TR) = (D{e). (3.1) 

For the purposes of this section, we do not need to make other assumptions. The gauge algebra 
may be irreducible or reducible, and close off shell or on shell. The theory may be renormalizable or 
nonrenormalizable, perturbatively unitary or not. In particular, it may contain higher-derivative 
fields. The action S does not need to satisfy special cohomological properties. We can also include 
local composite fields 0^ (x), in renormalizable and nonrenormalizable theories, by coupling them 
to external sources Lj{x) and appropriately extending the actions Sc, Sd, Sd and S. In the 
arguments that follow, the dependence on such types of external sources is not made explicit. 
However, we understand that it may be there, whenever necessary. 

Consider a canonical transformation d>, AT ^ ^',K' with generating functional 

Fi<^,K',e) = j + (3.2) 

where Q = 0{6) is a local functional. Let Sg denote the action obtained by applying (3.2) to S, 
Sr$ the renormalized version of S$ and TRg the renormalized F functional associated with Srq. 
We assume, for simplicity, that Q does not contain analytically evanescent contributions. 

We work out how F^ and the identity (3.1) change when we make the transformation (3.2) on 
S. To reach Sg from S, it is useful to embed the theory into a more general theory, by considering 
the extended action 

= S{^,K,u) + (3.3) 

i 

where r* are arbitrary parameters and {Hfj} is a basis of local functionals of <I> and K. Specifi¬ 
cally, the dCj are integrals of local monomials constructed with the fields, the sources and their 
derivatives. They can be restricted by demanding that they be invariant under the nonanoma¬ 
lous symmetries of the theory. However, they are not restricted by gauge invariance, or power 
counting. To simplify a number of formulas, we include duplicates of the terms that are already 
present in S, multiplied by new independent parameters tiTi. The difference S — S' is made of 
0(^)-terms and is also assumed to contain evanescent terms (including those that are already 
present in S). Basically, S — S parametrizes the arbitrariness of the subtraction scheme. We 
denote the F functional calculated with the action S by n(<h, AT, w, ^r). 

Now, we renormalize S. We denote its renormalized action by Tir and the F functional 
associated with Tir by Q,r. We can imagine, for a moment, that we replace each Kti with an 
ordinary parameter pi of order zero in h. In that case, the construction of S/j is straightforward. 
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since every divergence can be subtracted by means of pi redefinitions. At a second stage, we raise 
the order of the parameters pi by restoring in their places. The consistency of this operation 
is justified by the arguments that follow. 

We organize the renormalization of S so that coincides with Sr when the parameters r* 
are equal to suitable finite functions which identify the subtraction scheme where formula 

(3.1) holds: 


j:R{^,K,u,hT*) = SR{^,K,u), nR{^,K,u,hT*) = rR{^,K,u). (3.4) 

At arbitrary r, the action T,r can be viewed as an extended renormalization of S, which includes 
the most general subtraction scheme. We say that 'Sr is the arbitrary renormalization of S. When 
we set Ti = T* we specialize the subtraction scheme to the one used for Sr, which, by assumption 
(3.1), preserves gauge invariance to all orders. 

Since it is consistent to set r* = r*, it is also consistent to set Ti = t* + K^i'n+ii, n ^ 0, 
for arbitrary new parameters hn+ii- By this we mean that the renormalization of each hn+ii 
remains analytic in h. We can better explain this fact by noting that the renormalizations of the 
differences 6i = r* — r* vanish at 6j = 0, so they must be proportional to 5j. Thus, if we replace 5i 
by h^Hn+ii, n > 1, the renormalizations of izn+n remain analytic in h. These remarks illustrate a 
trick that we use in the recursive proof given below. Precisely, at each step we raise the h order of 
certain residual parameters by one unit, till we make those parameters disappear, and show that 
we can do this while preserving the analyticity in h. 

The dehnition (3.3) understands that the difference S — S starts from 0(H). Indeed, we do 
not want to modify the classical action, but just parametrize the arbitrariness of the subtraction 
scheme. The reason why we move to the more general theory S is that if we want to cancel the 
anomalies after the canonical transformation, we generically need to re-hne-tune all sorts of finite, 
local terms, including the gauge noninvariant ones. 

As said, S$(^, K,ui,9) denotes the action obtained by applying (3.2) to 5'(<h, iP, cj). Let 
Sg(^, K,u},hT,6) denote the action obtained by applying (3.2) to S. We obviously have Sg = 
Sg + 0(h). We denote the renormalized version of Sg by SRg(^,K,uj,hT,9). Since SRg = 
Sg + 0(H) = Sg + 0(h), SRg Can be viewed as the arbitrary renormalization of Sg. Note that Sg is 
not gauge invariant, so its renormalization is not subject to particular restrictions, aside from the 
continuity condition SRg(^,K,oj,hT,9) = Sr(^, K,uj,hT) + 0(0). We denote the T functional 
associated with SRg by il.Rg. 

Finally, consider the local functional Q(^,K') defined by the canonical transformation (3.2), 
and dehne Qg(^,K') = dQ(^, K')/d9 and Qg(^,K) = Qg(^, K'(^, K)). Let QRg denote the 
renormalized version of Qg(^,K) at generic r. 

We prove that 
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Theorem 1 there exist finite functions = 0{9), such that, defining f*(uj,9) = Tj{^) + 

t'*{ijJ,9), the action 

SRe{^,K,co,9) = J:Re{^,K,uj,hf*{uj,9),9) (3.5) 

gives a T functional that satisfies the identities 

= 0(e), (3.6) 

- (r^e, (<5/?e)) = 0(e), (3.7) 

for arbitrary 9, where Qrq denotes the functional Qrq calculated at hxi = hff. 


Note that formula (3.5) ensures that Srb also satisfies the continuity condition SRe(<h, K, ui, 9) = 
Sr{^, K,uj) + 0{9). In fact, all the operations we make preserve the continuity in 9. 

For clarity, it is useful to summarize the dehnitions given so far in a table: 


S 







Sr 

r 

—)• 


fhr* 


fhT* 

^R 

r 

— 

^R 




'S‘RjB 

r 

- 

^Re 




Srb 

r 

- 

r_R6» 


3.1 The equations of gauge dependence 

If we apply the identity (A.5) of appendix A to the renormalized action Tirq and the renormalized 
r functional with X = Qro, we obtain 

- {^Re, {Qro)) = ^ - (S_R6», Qrb) - ^C^Rd, S/J6))Qr6)^ . (3.8) 

It is convenient to organize this formula in the form 

= {O,R0, (Qro)) + (^ii;6>)r > (3-9) 

where 

^Rd = --^i^Re, ^Re)QRe H —- {'^Re, Qro)- (3.10) 

If the right-hand side of formula (3.9) contained no {^Re)^ (which happens, for example, in 
manifestly nonanomalous theories) or we knew that {^rb)^ is for some reason equal to 0(e), the 
solution of our problem would be straightforward. Formula (3.9) would turn into a much simpler 
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equation, which is integrated in ref. [ 8 ] and in appendix C. The result would be that the entire 
6 dependence of fl/jg can be absorbed into a convergent canonical transformation acting on 
up to 0(e). Moreover, there would be no reason to keep r generic. More simply, we could just 
work with r = r* from the start. Then, formula (3.9) would give (3.7). Integrating (3.7) with 
the procedure of appendix C, we would find a convergent canonical transformation that turns T^^ 
into T/je, again up to 0(e). That canonical transformation would also turn formula (3.1) directly 
into (3.6), since the right-hand side would remain evanescent. 

Unfortunately, (V_R 6 »)r i® there, because the theory we are considering is potentially anomalous, 
so we must study the effects of such an extra term. To achieve this goal, a few facts need to be 
noticed. 

(i) By construction, Qjig and (Qro) are convergent. 

(a) The local functional {'ERg,T,Rg) is already renormalized. Indeed, formula (2.9) tells us 
that ((S/je, S/je)) = (Il/je, fl/je), which is convergent. Since = Sg + we can say that 

’^Rg) is the arbitrary renormalization of {Sg,Sg). 

(in) By points (i) and (ii), all the subdiagrams of the diagrams that contribute to the av¬ 
erage {{T,Rg,'ERg)QRg)Y are already renormalized, except those that contain both insertions of 
(S_R 0 , T,Rg) and QRg. 

{iv) The object ^Rg is a bit peculiar, because at the tree level it is equal to 

Yg = -'-{Sg,Sg)Qg. (3.11) 

The reason why the last two terms of (3.10) do not contribute at Ii = 0 is that 

= = 0, (3.12) 

which follows from formula (A. 6 ), if we understand that the primes denote the fields and the 
sources before the transformation, i.e. write S = and Sg = Sg{^,K). We see that Yg is 

the product of two local functionals. We call Yg a local bifunctional. We extend the definition of 
local bifunctional to any expression of the form 

T, = '^AiBi + C (3.13) 

i 

where Aj, Bi and C are local functionals. An evanescent local bifunctional is a local bifunctional 
(3.13) where C and Aj (or Bi) are evanescent. 

Now, {S,S) is an evanescent local functional, by formula (2.7), and Sg is obtained from S by 
means of a finite canonical transformation, which preserves the antiparentheses and maps 0 (e) 
into 0(e). Thus, (Sg, Sg) is also evanescent, and ^Rg is an evanescent local bifunctional. Actually, 

(v) is a renormalized evanescent local bifunctional, since formula (3.9) implies that (y_R 0 )p 
is convergent. 
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The procedure to renormalize a local bifunctional is explained in appendix B. There, it is also 
shown how to renormalize an evanescent local bifunctional £ in such a way that (£_r)p = 0(e)- 
To describe what happens order by order in the perturbative expansion, consider for simplicity 
an evanescent local bifunctional of the form £ = EB + F where E and F are evanescent local 
functionals. Let En and Bn denote the functionals E and B renormalized up to and including n 
loops, and inductively assume that En satishes {En) = 0(e) + Also assume that Fn is a 

local functional such that (£n)r = 0(e) + where E-n = EnBn + Fn- Then, the 

contributions to (£n)p are the sum of a local divergent part, a local nonevanescent part and a 
generically nonlocal evanescent part. If Bn+i is the functional B renormalized up to and including 
n + 1 loops, there exist local functionals En+i and Fn+i such that (En+i) = 0(e) + 0(li"'+^) and 
(£n+i)p = 0(e) + where £n+i = En+iBn+i + Fn+i- The subtraction can be iterated in 

n to obtain (Er) = 0(e) and (£_r)p = 0(e), where Er = E^o and Er = £oo- 

Although is renormalized, it does not satisfy (y_R6»)p = 0(e), as far as we know. However, 
we will obtain = 0(e) by identifying the functions 9) and setting r* = t* + 

To prove (3.5), (3.6) and (3.7), we proceed by induction. Let Vnj denote free parameters of 
order fE. The hrst inductive assumption is that 

(a„) there exist finite functions |U„j(ti;, zz^+ifc, ^) = 0{9)0{h), such that the action 

Sn(‘h, iL, CU, Z^n+lj ) ^) = (‘h, iL, CU, h-Tj finj E l^n+lj ^ 9) (3.14) 

gives a T functional that satishes 

{Qn, ^n) = ((Sn, S„))„ = 0(e) + 0(h"+^), (3.15) 

where (•••)„ denotes the average calculated with the action 
Now, dehne 

Qn — Qr9 (*^) ^ ■) j T hni T ^n+lj j9), 

i dT, 

Yn = --{^n,^n)Qn + ^-{^n,Qn)- (3.16) 

Applying formula (A.5) to the action and its T functional with X = Qn, we obtain 

QQ — (^n; {Qn)n) + (^n)„p , (3-17) 

where (■ ■ ■ )nr denotes the one-particle irreducible diagrams of the average (■ ■ • )n- second 

inductive assumption is that 

(bn) 

(F„)„P = 0(e) + 0(r+i). (3.18) 

Statement (oq) is true with = 0, because Sq = S'g -|- 0(h) and {So,S$) is evanescent, so 
((So,So))o = 0(e)-|-0(h). Statement (6o) is also true, because Qo = Qe + 0(h), so Lq = h'6» + 0(h). 
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3.2 Inductive proof 


Assume that (an) and (bn) hold. Then, the averages ((S„, 'En))n ^i^cl (Xn)nT evanescent up to 
and including n loops. The arguments of appendix B ensure that the (n + l)-loop contributions 
to (Xn)nT^ which are convergent by formula (3.17), are the sum of a local nonevanescent 
part Y^onev plus a generically nonlocal evanescent part. We have 


(y„)„r = 0 (e) + yii+ii + 0(r+2). 


(3.19) 


We can write an explicit expression for yln(iiit. Recall, from formula (3.3), that the derivatives 
dY/d{hTj) form a basis for the local functionals of and K. Obviously, so do the derivatives 
dYg/d{hTj) = “Kjo. Up to higher orders in h, the derivatives dYjig/d{hTj) = “KjQ + 0{K) are also 
a basis, as well as the derivatives dYn/dvn+ij- Thus, there exist finite order-^'^'*'^ functions 
which depend analytically on ca, I'n+ik and 6, such that 


y(n+l) _ 
nnonev / ^ 


(J 


{n) dYn 


^ dVn+lj 


+ 0(^’^+^). 


Now, define 


^n+i = yn - O' 


(n) 


dYr, 


^ dVn+lj 


Taking the average of both sides, and using (A.3), we get 


(3.20) 


(3.21) 




i 




dv. 


11+1 j 


(3.22) 


Using (3.19) and (3.20), we obtain 

{^n+l)nT = ytinl ” Yttnl + 0(e) + 0(^+2) = 0{e) + 0(^+2). (3.23) 


Using (3.22) inside (3.17), we also find 

^ = (!!„, ((3„)„) + y + <«»+i>nr ■ (3-2'l) 

Define finite functions i/„_|_ij(a;, 0) as the solutions of the evolution equations 

=-<jf\uj,un+ik,e), (3.25) 

with the initial conditions r'„_|_ij(a;, Pn+ifc; 0) = Un+ij- Clearly, Vn+ij = i^n+ij + 0(0). Given a 
functional A(d>, it', w, 0), define 

A($, A',a;,P„+ij,0) = X(d>, AT, w, r'n+y (w, P„+ifc, 0), 0). (3.26) 
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Then 


dX 


dX _(„) dX 
do ^ dun+ij ’ 


(3.27) 


where are the functions obtained by applying the redefinitions Vn+iji'^i^n+ik-,0) to 
Choosing X = we can turn equation (3.24) into 


^ = (^n, (gn)„) + (yn+l)„r. (3-28) 

Applying the redefinitions Un+ij{oJ,Vn+ik^O) to the functions ^„j(a;, r'yi+ifc, ^) of assumption 
(a,i), and including the contributions coming from Vn+ij — which are proportional to 6, we 

can define new 0{0)0{h) functions Jinj{oJ,i>n+ikiO) by the formula 


l-^nj (^) ^^n+1 fc) 0) — Unj (w, r'jj+l fc (^j ^n+li; ^) > ^) T ^n+lj i k'n+lkj 0) . 

Then, using (3.26) and (3.14), we have 

^n(^) at, U!, l^n+lj j O') — S/J0(4?, AT, CU, flTj + Unj + , O'). 

At this point, the independent parameters are w, T'n+ij and 0. The formulas we have written 
so far hold for every value of r'n+ij; as long as it is Now we want to raise the h order of 

i^n+ij by one unit. The validity of this choice will be self-evident. By this we mean that it allows 
us to iterate all the arguments of the proof without difficulties till the very end and preserve the 
analyticity in h. 

Define 


k'n+lj — k'n-\-2ji An+lj (^) 2^n+2fc) ^) — Anj (a^j 2^n+lfc; ^) • (3.29) 

So doing, we obtain the action Sjj+i, given by formula (3.14) with the replacement n ^ n + 1: 

^n+l (^) at, CU, l'n+2j i0) — ^Rj9 (*^) A", UJ , j -\- ^n+lj T ^n-|-2j }0') — Sjj ($, K, UJ , i0'). (3.30) 

Recalling that = S/j -|- 0(0) and fin+ij = O(0)O(^), formula (3.30) tells us that, at 0 = 0, 

Sn+l|,=o = ^R{^,K,U,hT; + l^n+2j) = , K, CO, flT*) + 0(^+2) = SR{^,K,io) + 0(^+2), 

where the last equality follows from the first equation of (3.4). Finally, the second equation of 
(3.4) and formula (3.1) give 

(D„+i,D„+i)|,^o = (rR,r^) + 0(r+2) = (D{e) + 0(r+2). (3.31) 

This is a check that the new action S„+i is AB nonanomalous at 0 = 0, up to 0(e) and 0(/i"'+^). 
Now we show that satisfies the same property for every 0. 
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Using formula (3.28), we get 


de 


(^n+l) {Qn+l)n+l) + (^+l)jj+ir ’ 


(3.32) 


where the functionals Qn+i and are obtained from Qn and ^n+i by applying the redefinitions 
Un+ij{oj, Un+ik, ^) and (3.29). Using (3.26), (3.27) and (3.21), it is easy to see that formulas (3.16) 
hold with n —)■ n + 1. 

Moreover, formula (3.23) ensures that 


(y„+i)„+ir = 0(e) + 0(r+2), 


that is to say formulas (3.17) and (3.18) hold with n —)• n + 1. 

Taking the antiparentheses of (3.32) with fln+i and using the Jacobi identity, we also find 

d 

(f^n+l) f^n+l) — ((f^n+1) f^n+l)) (Qn+1 )n+l) T ^(nn-i-i, (yj-|_l)^_l_j^p). (3.33) 

The last term of the right-hand side is 0(e) + In appendix C we show how to integrate 

equation (3.33) and prove that the 0 dependence of (fln+ijUn+i) is encoded into a canonical 
transformation, up to 0(e) and By formula (3.31), the value of (fln-i-i; ^n-i-i) at 0 = 0 is 

also of such orders. Moreover, the canonical transformation is convergent, because it is uniquely 
determined by {Qnj^i)n+i, which is convergent. Therefore, we find 

(fl„+i,0„+i) = 0(e) + 0(r+2) 


for arbitrary 9, which is formula (3.15) with n —n -|- 1- As promised, the action S„_|_i is AB 
nonanomalous for arbitrary 0, up to 0(e) and We have thus proved statements (a^+i) 

and {hn+i)- 

Finally, formulas (3.6) and (3.7) follow by taking n to infinity, with Vooj = 0 and hTj'{u!,9) = 
f^ooj{oj,0,6). Indeed, because of (3.14), if we define Sj^ according to (3.5), we have Sqo = Sjio, 
so rioo = T/jg). Then, formula (3.15) becomes (3.6) at n = oo. By (3.18), formula (3.17) turns 
into (3.7) at n = oo, with Qjiq = 


3.3 Integrating the equations of gauge dependence 

Equation (3.7) can be integrated with the method of appendix C (see also [8]). There, it is shown 
that we can consistently ignore the terms 0(e) appearing on the right-hand side, in the sense that 
the solution we find by ignoring those terms is correct up to 0(e). The basic reason is that the 
equations involve only convergent functionals. Alternatively, we can just remove the cutoff by 
taking the physical limit e —)■ 0 in (3.7) and then work in the physical dimension d. The result is 


21 



that every 9 dependence of F/jg can be absorbed into a convergent canonical transformation, up 
to 0(e). 

More precisely, the theorem of appendix C ensures that there exists a canonical transformation 
<f>, iF —)■ <h', iF' such that the F functional F^ dehned by 

F'^(d>', K', u) = F^,(d>(<f>', K', uj, 0), iF(d>', iF', a;, 0), cu, 9) (3.34) 

is 9 independent, up to 0(e). Setting 0 = 0, we hnd F):j = F^, since 

F^(ch', iF', uj) = F^0(<h', K\ w, 0) = F^(d>', K', u). 

Finally, inverting the transformations, we get 

Fnei^, iF, w, 9) = FR(d>'(d>, K, co, 9),K'{^, iF, w, 0), w). (3.35) 

As promised, the dependence of F/j^ on the gauge parameter 9 can be fully absorbed inside a 
canonical transformation. 

We recall that the canonical transformations we are talking about, which are convergent, 
nonlocal and act on the renormalized F functional, originate from a local canonical transformation 
of the form (3.2) that acts on the tree-level action. The connection between the two is a procedure 
of re-renormalization and a re-hne-tuning of the hnite local counterterms. We call such canonical 
transformations on F^j special. Clearly, the composition of special canonical transformations is a 
special canonical transformation. If we repeat the argument of this subsection for any other gauge 
parameter 9 that satishes (3.7), taking one at a time, we can prove that the entire dependence of 
the F functional on the gauge parameters can be absorbed into a special canonical transformation. 

In subsection 7.3 the equations of gauge dependence are used to prove that the physical 
quantities are gauge independent. 

4 Gauge dependence of the renormalized action 

In this section we study the counterparts of equations (3.6) and (3.7) at the level of the renor¬ 
malized action. Using the identity (2.9), formula (3.6) gives (Fij6i,F^6i) = {{Sne, Sro)) = 0(e), 
which implies that {Srq,Srq) is a “truly evanescent” local functional, i.e. a local functional such 
that its average is evanescent. We use the symbol £ to denote such type of functionals. Thus, we 
have the formula 

iSRe,SRe) = ^, (4.1) 

where (£) = 0(e). Equation (4.1) expresses the cancellation of the gauge anomalies to all orders 
at the level of the renormalized action. 
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Next, if we apply formulas (A.3) and (A.4) to (3.7), we obtain 

^ - -:^{{SRe-,SRe)QR 0 )Y = C)(e). (4.2) 

By formula (4.1), {Srq, Srq) is a renormalized local functional such that its average is evanescent. 
In appendix B we prove that there exists an 0(h) local functional Fr, such that the local bifunc¬ 
tional Yji = —{i/‘i){SR 0 , Sr 0 )Qr 0 + Fji is renormalized and the average is evanescent to 

all orders. We denote such Fr by the symbolic expression {Sr 0 , Sr 0 ;Qr 0 ). Thus, formula (4.2) 
gives 

^ - iSR 0 ,QR 0 ) - iSR 0 ,SR 0 -,QR 0 )^ = 0 (e). 

In turn, this equation implies 

dSR 0 

QQ = {Sr 0 , Qro) + {Srjb, Sr 0 ] Qr 0 ) + £. (4.3) 

Formula (4.3) is the equation of gauge dependence for the renormalized action Srq. Note that 
{Sr 0 ^ Sro; Qr 0 ) encodes the re-hne-tuning of the finite local counterterms. 

Equation (4.3) can be integrated with the method explained in appendix C. Although the term 
{Sr 0 , Sr 0 ; Qr 0 ) depends on Sr 0 , a recursive procedure allows us to treat it as a known functional 
at every step. 

5 RG invariance and other applications 

In this section we give a few applications of the theorem proved in section 3. The first application is 
the proof that RG invariance is preserved by the canonical transformation. The second application 
is the proof that renormalizable chiral gauge theories gauge-fixed by means of a nonrenormalizable 
gauge hxing remain renormalizable, although in a nonmanifest way. The third application is a 
step of the proof of the Adler-Bardeen theorem in nonrenormalizable theories [7]. 

RG invariance is expressed by the Gallan-Symanzik equation (which is derived at the end of 
this section) 

where /3* are the oji beta functions (at e 7 ^ 0) and Ur is a local functional. At the level of the 
renormalized action Sr, the Gallan-Symanzik equation reads 

(^- 2 ) 

a/i ouji 

Let fij denote the renormalized T functional where the parameters uii are written in terms of 
their running versions and where tD* are the solutions of ^idoji/d^ = with initial 


23 



conditions Wj. We have 


Let Ur denote the functional (Ur) reparametrized in a similar way. Then the Callan-Symanzik 
equation becomes 

^^-{rR,UR) = 0 {e), 

OfJ, 

where d denotes the derivative at fixed A. The new equation has the same form as (3.7), so it is 
solved by making a canonical transformation. From formula (3.34), we learn that there exists a 
canonical transformation that takes us to new fields and sources <L, K and a reference value of fj,, 
which we denote by jl and leave implicit, such that 

for a certain other functional T/j. 

Now, if we make the canonical transformation (3.2) on the tree-level action, we get, by formula 
(3.35), 

f^0(<h, K, oj, fi, 6) = fjj(d>'(4>, K, Cj, li, K, uj, fi, 0), w, m) = 

= f H(4>(d>', iL', cD, m), ^(<h', iL', w, m), w). 

Going back to the parameters u, we also have 

VRei^, K, u, fi, 9) = fR,(<L, K, Cj, 9) = f ^(l>'(<h, K, uj, 9),K'{^, K, u, 9),oj{io, fi)). 

having defined <L(<h', K', oj, n) = ^'(<h, K, oj, 9) and similarly for K. Differentiating with respect 
to In fi, and recalling that our canonical transformations are special, we get 

- (Xro, (URe)) = 0(e), (5.3) 

for some new local functional Urq. Formula (5.3) is the transformed RG equation. Note that the 
beta functions do not depend on 9 in this approach. 

Another application that we mention is to power-counting renormalizable chiral gauge theories 
gauge-fixed by means of a nonrenormalizable gauge fixing. If a renormalizable theory is nonchiral, 
it is rather straightforward to prove that it remains renormalizable when a nonrenormalizable 
gauge fixing is used. When the theory is chiral, on the other hand, the matter is more complicated. 
In principle, the simplifications between divergences and evanescences can make the parameters 
of negative dimensions, introduced by the gauge fixing, propagate into the physical sector and 
turn the theory into a truly nonrenormalizable one. The theorem of section 3, combined with RG 
invariance, ensures that this cannot happen. 
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Consider for example the standard model in flat space and gauge fix a non-Abelian gauge 
symmetry by means of a gauge-fixing function such as 

Since the constant k has dimension —2 in units of mass, power counting alone is not sufficient 
to classify the counterterms in a convenient way at k 7 ^ 0. However, the change of gauge fixing 
that turns G“( 0 ) = d^A^ into G“((/)) is a canonical transformation, so we can apply the theorem 
of section 3. Formula (3.3) teaches us that infinitely many terms IKj of arbitrary dimensions are 
switched on, including the gauge noninvariant ones. Nevertheless, the theorem ensures that once 
we have done that, it is possible to express the coefficients of all of those terms as functions of the 
other parameters of the theory, and fine-tune those functions to enforce again the cancellation of 
gauge anomalies to all orders. Moreover, the argument given above ensures that RG invariance 
is preserved. We conclude that no new independent parameters are necessary to subtract the 
divergences and cancel the gauge anomalies in a RG invariant way: the physical couplings are 
still finitely many. Thus, when a power counting renormalizable chiral gauge theory, such as 
the standard model in flat space, is gauge-fixed by means of a nonrenormalizable gauge fixing, 
it remains a renormalizable theory, although its renormalizability is not manifest anymore. A 
similar conclusion holds when the theory is renormalizable by weighted power counting [27] or 
any other criterion. 

The third application we mention is the proof of the Adler-Bardeen theorem in nonrenormal¬ 
izable theories, recently obtained in ref. [7] by upgrading the arguments of [ 6 ]. It applies to the 
theories whose gauge symmetries are general covariance, local Lorentz symmetry and Abelian and 
non-Abelian Yang-Mills symmetries, and satisfy a variant of the Kluberg-Stern-Zuber conjecture. 
Quantum gravity coupled to the standard model satisfies all the assumptions and so is free of 
gauge anomalies to all orders. In the approach of [7], the GD regularization is combined with 
a higher-derivative regularization. If the scale A associated with the higher-derivative terms is 
kept fixed, we obtain a super-renormalizable higher-derivative (HD) theory, which satisfies the 
Adler-Bardeen theorem by simple power-counting arguments. When the scale A is sent to infinity, 
the A divergences are renormalized inductively. At each step, the theorem of section 3 allows us 
to resubtract the divergences in e and re-fine-tune the finite local terms, in order to enforce the 
cancellation of gauge anomalies to all orders at A fixed. In the end, thanks to this, the cancel¬ 
lation of gauge anomalies survives the renormalization of both types of divergences. Moreover, 
the approach of ref. [7] identifies a special subtraction scheme where the cancellation of gauge 
anomalies is manifest from two loops onwards, within any given truncation. We stress that it is 
not possible to achieve a similar goal by means of regularization-independent methods. 

The Gallan-Symanzik equation (5.1) can be proved from the results of ref. [7] as follows, 
under the assumptions specified there. At A fixed the HD theory is renormalized by redefinitions 
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of parameters, while the trivially anomalous terms are canceled by adding a finite local functional 
—x/2 to the action. The renormalized action coincides with its bare version apart from y itself, 
which satisfies y = where XB is RG invariant. Then, the HD theory satisfies formulas 

(5.1) and (5.2) with Ur = 0, the right-hand side of (5.2) being equal Exl‘^- The average (y) 
is convergent in the HD theory (since its divergences, which would start from two loops, are 
excluded by the arguments of [7]), so the product e(y)/2 is truly evanescent at A hxed. At a 
second stage, the renormalization is completed by removing the A divergences. This is done by 
means of special canonical transformations and redefinitions of parameters. In this section we 
have proved that those operations preserve the Callan-Symanzik equation, although they can 
affect the beta functions and the functional Ur. In the end, we obtain equations of the forms 
(5.1) and (5.2). 

6 Gauge dependence of the beta functions 

Often, we can prove that a theory is AB nonanomalous in a family of gauges, parametrized by 
certain gauge-fixing parameters In various common situations we can achieve this goal by 
applying the results of ref. [6], where the Adler-Bardeen theorem was proved for arbitrary values 
of the gauge-fixing parameter ^ of the Lorenz gauge, in power counting renormalizable gauge 
theories that have unitary free-held limits. More generally, if the theory is coupled to quantum 
gravity, we can apply the results of [7]. Then, when we study the dependence of the correlation 
functions on we can proceed more straightforwardly than in section 3, since we already know 
that (rR,r^) = 0(e) for arbitrary It is worth recalling that in section 3 we had to derive this 
result from just knowing that (r/j,r/j) was 0(e) for ^ equal to some initial value 

In this section we study the equations of gauge dependence in theories that are AB nonanoma¬ 
lous for arbitrary values of some gauge parameter 6 and satisfy some additional assumptions. 
Those assumptions are not very restrictive, since they are fulhlled quite commonly. When 9 
varies, we do not need to readjust the subtraction scheme by hne-tuning the hnite local countert¬ 
erms. Then, however, the beta functions of the couplings are in general gauge dependent. Their 
gauge dependence can be removed by redefining the couplings themselves. 

We begin by listing the assumptions we need. 

(I) We assume that the gauge algebra is irreducible and closes off shell. This assumption is 
satisfied by the theories whose gauge symmetries are general covariance, local Lorentz symmetry 
and Abelian and non-Abelian Yang-Mills symmetries, such as the standard model coupled to 
quantum gravity. It allows us to make a number of simplihcations. For example, we can choose 
the fields and the sources K so that the gauge-fixed tree-level solution Sd of the D-dimensional 
master equation (Sd, Sd) = 0 is linear in K and has the very simple structure (2.3). 

We have already remarked that in various cases, for example when the theory is chiral or parity 
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violating, the action Sd, embedded in D dimensions using the standard rules of the dimensional 
regularization technique, is in general not well regularized, due to the key role played by the 
d-dimensional analogue 7 of the matrix 75 , or the tensor Using the chiral dimensional 

regularization, a well-regularized classical action S{^,K) is obtained by adding a number of 
evanescent corrections S'ev to Sd [11], as shown in formula (2.6). We denote the parameters 
contained in S'ev by r]!- For convenience, we assume that S'ev depends linearly on the parameters 
r], and vanishes for 77 = 0 . 

Let {Si(</')} denote a basis of local gauge invariant functionals of the classical fields (j). Expand 
the classical action as 

5e(</.) = ^A,g,(</.), (6.1) 

i 

where A, are independent parameters. We call the constants Aj “physical parameters”, since they 
contain or are related to the gauge coupling constants, the masses, the Yukawa couplings, etc. 
In our notation some parameters A* may be actually redundant. Nevertheless, to simplify some 
derivations we prefer to keep an independent Aj for every Sj. For example, it is often useful to 
restrict Sc by dropping the terms that are proportional to the Sc field equations, because those 
terms can be renormalized by means of canonical transformations, rather than A* redefinitions. 
We do not implement this restriction right now, to make some arguments of the derivations that 
follow more transparent. We can always remove that class of redundant terms at the end by 
means of a convergent canonical transformation, by applying either the procedure of section 3, 
which is more general, or the one of this section, which holds under specific assumptions. Both 
procedures preserve the cancellation of gauge anomalies and the equations of gauge dependence. 

In total, we have physical parameters A, gauge-fixing parameters contained in and 
regularizing parameters r/. The classical action is written as S{^,K,X,^,r]). 

The action Sc may contain accidental symmetries, which are the global symmetries unrelated 
to the gauge transformations. Some accidental symmetries are dynamically lost, because they 
are anomalous, others are nonanomalous. Let G^as denote the group of nonanomalous accidental 
symmetries, or the identity group, depending on whether the gauge group contains U{1) factors 
or not. By definition, the set {9i{(p)} includes the invariants that explicitly break the anomalous 
accidental symmetries, but excludes the invariants, denoted by Si(</>), that explicitly break Gnas- 
Then the actions Sc and Sd do not contain the invariants Sj, so we define extended actions Sc 
and Sd = Sc + {Sk,'^) + Sk that do include them, multiplied by independent parameters A*. 
Both choices of including and excluding the invariants Sj, are consistent, from the point of view 
of renormalization. 

We say that the action Sd satisfies the Kluberg-Stern-Zuber assumption [19], if every nonevanes- 
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cent local functional X of ghost number zero that solves the equation {Sd,X) = 0 has the form 

X = ^ai9i + {Sd,Y), 

i 

where a* are constants depending on the parameters of the theory, and y is a local functional of 
ghost number —1. 

We say that the action Sd is cohomologically complete if its extension Sd satisfies the extended 
Kluberg-Stern-Zuber assumption, that is to say every nonevanescent local functional X of ghost 
number zero that solves {Sd,X) = 0 has the form 

X = Y,ai9i + Y,bi9i + {Sd,Y), (6.2) 

i i 

where bi are other constants, and y is a local functional. 

(II) We assume that the action Sd of (2.3) is cohomologically complete and the group Gnas is 
compact. 

The Kluberg-Stern-Zuber assumption is satisfied when the Yang-Mills gauge group is semisim¬ 
ple and the action Sd satisfies generic properties [30]. It is not satisfied when the gauge group 
has U{1) factors and accidental symmetries are present. In particular, it is not satisfied by the 
standard model. However, it can be proved, using the Ward identities that hold in the Lorenz 
gauge, that the standard model is cohomologically complete [6]. So are the Lorentz violating 
extensions of the standard model of refs. [12, 31], which are renormalizable by weighted power 
counting [27]. Starting from the cohomological theorems proved in ref. [30], it can be proved that 
the standard model coupled to quantum gravity is also cohomologically complete [7], and so are 
most of its extensions. 

The condition (Sd, X) = 0 is the one typically satisfied by the counterterms. In this section we 
show that the contributions of the extra term contained in the generalized Ward identity (2.15) 
satisfy the same condition. Thus, assumption (II) will give us control on the effects of the new 
term. 

We can imagine that 9 is one of the parameters or another parameter introduced by a field 
redefinition. We keep it distinct from the other parameters X,^,r] contained in the action S and 
assume that S{^,K,X,^,r]) denotes the action at some specific value 9* of 9. With no loss of 
generality, we take 9* = 0. By definition of gauge parameter, when we vary 9, we make a canonical 
transformation generated by a functional of the form (3.2) on the action S, and this operation 
gives the action Sg. As before, let Sr$ denote the renormalized action and T/jg the T functional 
associated with it. 

(III) We assume that the theory is AB nonanomalous for arbitrary values of some gauge 
parameter 9. Precisely, we assume that there exists a class of subtraction schemes where the 
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renormalized T functional F/jg satisfies the identity 

(r_R 6 i,rij 0 ) = 0(e), (6.3) 

where 6 takes values in some continuous range that includes 9 = 0. From now on we understand 
that we work in that class of subtraction schemes. 

Assumption (III) has been proved, for common families of gauge conditions, in the power¬ 
counting renormalizable gauge theories that have unitary free-field limits [ 6 ], in the Lorentz vi¬ 
olating extensions of the standard model that are renormalizable by weighted power counting 
[12, 31], in the standard model coupled to quantum gravity and a large class of other nonrenor- 
malizable theories [7]. 

We prove that there exist finite functions pj of A, rj and 9, which start from 0{h), and a 
renormalized local functional = Q 0 {^, K) where Qg{^,K) = Qg{^, K'{^, K)) 

and Q 0 = dQjd9, such that F/jg satisfies the equation 

—^ = '^Pj-Q^ + i^Re,{HR9)) + 0{s). (6.4) 

j ^ 

The first term on the right-hand side of (6.4) can be absorbed by means of finite redefinitions 
of the parameters A (which correspond to the re-fine-tuning of the previous section). The second 
term is the one that can be absorbed into a canonical transformation. 

In the rest of this section we derive the equations of gauge dependence (6.4) under the as¬ 
sumptions listed above, and integrate them. Before beginning the derivation, a few preliminary 
remarks are in order. If we differentiate (6.3) with respect to any parameter we find 

=0(e). (6.5) 

Now we take the antiparentheses of both sides of formula (3.7) or (6.4) with F/jg, and use (6.5) 
for C = ^ s-iid C = the Jacobi identity satisfied by the antiparentheses and formula (6.3) again. 
At the end, we find a consistent relation of the form 0(e) = 0(e). Thus, we can view formulas 
(3.7) and (6.4) as the solutions to the condition (6.5) for C = 9. 

To explain this issue more clearly, let us define an operator Jr that acts on a (generically 
nonlocal) functional Y by taking its antiparentheses with F/j^: 5rY = (F/je,y). Formula (6.3) 
ensures that Jr is nilpotent up to 0(e), because the Jacobi identity gives 

J 2 y = (FKe,(F^ 0 ,y)) = ^i{TR 0 ,Tii 0 ),Y) = 0(e). (6.6) 

Therefore, it is meaningful to study the cohomology of Jr. Consider the problem JrT = 0, of 
which the e ^ 0 limit of (6.5) is an example. It is a nonlocal upgrade of the more standard 
cohomological problem {Sd,X) = 0, where X is local. Formula (6.5) tells us that dTf{ 0 /d 9 is 
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closed, in the sense of the 5r cohomology, up to 0(e). On the other hand, formula (6.4) ensures 
that there exist finite linear combinations of dTug/dC, that are (5r-exact, up to 0(e). 

However, nonlocal cohomological problems are difficult to solve and must be treated with 
care, because if we do not specify which nonlocalities are allowed and which are not, any closed 
functional can in principle be exact. In other words, we cannot derive (6.4) immediately from 
(6.5), which is why gauge dependence deserves a separate investigation. 


6.1 The equations of gauge dependence 

We apply formula (A.5) of appendix A to the renormalized action Sro and the renormalized F 
functional F/jg, with X = Qrq^ where Qrq denotes the renormalized version of the functional 
Qg{^,K). We obtain 

^ (Qro)) + {Ur0)y , (6.7) 

where 

Uro = --{Sro, SRe)QRe ^- {Srq, Qro). ( 6 . 8 ) 

Taking the antiparentheses of both sides of (6.7) with F^^ and using (6.5) and ( 6 . 6 ), we obtain 


(ri?0,(t/Ke)r) = O(e). (6.9) 

Differently from (6.5), this nonlocal cohomological problem can be reduced to a local one, and 
solved. The reason is that Urq is originated by an evanescent local bifunctional. We prove 
that there exist finite functions pj = 0{h) of A, rj and 6, and a renormalized local functional 
Wro = 0(h), such that 

{Um)r = (6.10) 

j ^ 

We proceed by induction. Assume that there exist finite functions p^j = 0(h) of A, rj and 
6, and a renormalized local functional Wn = 0(h), such that the partially subtracted functional 

Un = Uro — ^ Pnj Qy — (Sro, Wn) — -(Srq, SRg)Wn, (6-11) 

j ^ 

satisfies 

(Un)r = 0(e) + 0(h^+^). ( 6 . 12 ) 

This assumption is clearly satisfied at the zeroth order, where poj = 0 and Wq = 0, because by 
formula (3.12) we have Urq = Yg + 0(h), where Yg is evanescent and given by (3.11). 

Using formulas (A.3) and (A.4), we obtain the average 

(Un)r = {URe)r - “ (^^e, (Wn)), (6.13) 

j ^ 
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which is clearly convergent. Consider the (n + l)-loop contributions to {Un)^- They are 

convergent, because so is the right-hand side of (6.13). Moreover, the inductive assumption (6.12) 
states that the average (17n)r evanescent up to and including n loops, while (6.3) ensures that 
{iSRe,Sji0)) = {Tr0, r^g) is evanescent to all orders. The arguments of appendix B ensure that 
the functional is the sum of a local nonevanescent part C/nnoniv plus a generically nonlocal 

evanescent part: 

{Un)r = uttnl + m + (6-14) 

Thus, using (6.13) and (6.14), we have 

/9T 

{URe)v = + 0(e) + 0(^+2). (6.15) 

j ^ 

Inserting this expression inside (6.9) and using (6.5), (6.6) and (6.3), we obtain 

{rRe,Uttnl) = 0{e) + 0{h'^+^). 

Taking the (n -|- l)-loop nonevanescent contributions to this formula, we find 

iSde,Uttnl) = 0, (6.16) 

where S^e is the action obtained by applying the canonical transformation (3.2) to S^- In deriv¬ 
ing the result (6.16), it is important to recall that the tree-level action (2.6) and the canonical 
transformation (3.2) do not contain analytically evanescent terms. In turn, Sq and Sde satisfy the 
same property, and S^e is the full nonevanescent part of Sg. 

Applying the inverse of the transformation (3.2) to equation (6.16) and letting denote 

the functional obtained from U^tnev, we get 

VlS+l) = 0. (6.17) 

At this point, we apply assumption (II). Let us imagine that instead of working with the 
classical action Sc we work with its extension Sc, which includes the invariants Sj that break the 
nonanomalous accidental symmetries belonging to the group Gnas- Similarly, we extend Sd to 
Sd, 5'ev to Sev and S = Sd + Sev to S. Every extended functional reduces to the nonextended 
one when we set A = r) = 0, where A and fj are the extra parameters of Sc and Sev, respectively. 
If we repeat the operations that lead to (6.17), we obtain an extended, nonevanescent local 
functional IJ^tnlv that satisfies {Sd,U^tnev) = 0. By assumption (II), the action Sd satisfies 
the extended Kluberg-Stern-Zuber assumption. Therefore, there exist finite order-Ii"'^^ constants 
^{n+i) ^ ^(n+i) ^ depending on the parameters, and a finite nonevanescent local functional 1/^"'"*'^^ 
of order such that 

uttniv = E + E + (Sd, 
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If we set A = fy = 0 in this equation, we obtain 

Uttnl = + (5,, 


(6.18) 


where are equal to and at A = ?) = 0. However, 

^nnonev and Sd are invariant under Gnas; while the functionals Sj are not. If we average on 
Gnas (which we can do, since Gnas is assumed to be compact), the Si disappear or give linear 
combinations of the invariants Si, and turns into some We obtain^ 


(jin+l) = cr 
^ nnonev / ^ ^ 


(n+l) dSd 

dXi 


+ {SdX 


(n+l)N 


(6.19) 


for some new constants We have used Si = dSd/dXj. At this point, we apply the canonical 

transformation (3.2) again, and note that, by formula (A.6) the difference between the transformed 
dSd/dXj and dSde/dXj is equal to {Sde,Xg) for some local functional Xg. In the end, we get 


/'/■(n+l) _ (n+l) lie 

^nnonev / ^ \ {^ddi Ig ) 


dXi 


for some new local functional of order Now, define 


( 6 . 20 ) 


Un+l — Ur0 — E Pn+lj Qy — {Sro, Wn+l) — -{Sro, SR0)Wn+l, 


Pn+lj — Pnj + \ Wn +1 — Wn + '^Rg'^ \ 

where are the renormalized versions of the functionals Using (6.11), we also have 


t/„+i = Un - - \{SRg,SRg)V^^-^^\ (6.21) 

j ^ 

Recall that S = Sd + Sev, which implies Sg = Sdg + 0(e) and Srq = Sde + 0(e) + 0{h). Taking 
the average of both sides of (6.21), and using (A.3), (A.4), (6.20) and then (6.14), we find 


{Un+i)r = {Un)r - uttnl + 0(e) + 0 (^+ 2 ) = 0 (e) + 0 (^+ 2 ), 


which extends the inductive assumption (6.12) to the order n+l. Formula (6.10) follows from 
formula (6.15) for n = oo, with pj = p^oj and WRg = Woo = 0{h). 

^If the terms proportional to the Sc field equations are dropped from Sc, the average on Gnas may generate 
them back. In the case of general covariance, local Lorentz symmetry and Yang-Mills symmetries, the average of 
Si may also affect besides the coefficients but the final result is still of the form (6.19). 
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Finally, using (6.10) inside (6.7), we get 


+ Wro)) + 0(e). 

j ^ 

This formula is equivalent to (6.4) with the identihcation Hrq = Qrq + Wrq. Observe that Hrq 
is another renormalized version of the functional Q 0 {^,K), and just differs from Qrq by a choice 
of subtraction scheme. 


6.2 Integrating the new equations and RG invariance 


Now we integrate the equations (6.4). We can easily absorb away the first term on the right-hand 
side by making finite redefinitions X{\', 6 ) of the parameters A. We choose functions \i{\',^,rj, 6 ) 
that solve the evolution equations 

dXi 


de 




with the initial conditions Xi{X '= A(. Using formulas (6.4) and (3.27), we obtain 

dTRe 


de 


= (bije, {Hro)) + 0(e), 


( 6 . 22 ) 


(6.23) 


where f/jg is related to Tro according to the definition X{X',^,r],9) [see (3.26) and the arguments 
given right after that formula]. 

Observe that equation (6.23) is equivalent to formula (3.7) of section 3. This means the 
redefinitions Xi{X',^,r],9) perform the re-fine-tuning of hnite local counterterms (automatically 
incorporated in the approach of section 3) that was missing so far in the approach of the present 
section. As in subsection 3.3, equation (6.23) can be integrated with the method of appendix C. 
We find that there exists a canonical transformation <h,iF ^ ^',K' such that the T functional 
defined by 

T'r{^', K', a', C, h) = K', A', C, h, 9),K{^',K\ a', e, r?, 9), A(A', t?, 0), V, 9) (6.24) 

is 9 independent, up to 0(e). Since 0 = 0 gives T^ = T/j, we also have 

r^(<I>', K', A', e, r?) = rH0(<h', K', A', C, h, 0) = ri?(<h', K', A', e, r?). 

Inverting the transformations, we obtain the formula 

TRei^, K, A, e, V, 9) = r^(<h'(<h, K, A, C, V, 9),K'{^, K, A, V, 9), A'(A, C, V, 9),^, r?), (6.25) 


which shows that the dependence of Trq on the gauge parameter 9 can be fully absorbed inside 
a hnite redehnition of the parameters A and a canonical transformation. 
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According to formulas (6.24) and (6.25), the beta functions I3'y of the parameters \' (in the 
framework where the helds and the sources have primes) are 9 independent. That means, however, 
that the beta functions I3\ of the couplings A do depend on 0. However, their 0 dependence is not 
arbitrary, because it disappears by making the redefinitions X{X',^,ri,9). 

We can repeat the argument for any other gauge parameter 6 for which formula (6.3) is known 
to hold, taking one at a time. Since the composition of special canonical transformations and re¬ 
definitions of parameters is a special canonical transformation combined with a redefinition of 
parameters, we reach the conclusion that the entire dependence on the gauge parameters can be 
absorbed into such operations, which do not affect the physical quantities (see subsection 7.3). 

We can also repeat the arguments of section 5 and prove that RG invariance is preserved. 
The difference is that now instead of (5.3) we get a transformed Callan-Symanzik equation that 
contains ^-dependent beta functions. 

7 Gauge independence and unitarity 

In general gauge theories we need to introduce extra fields, such as the Fedeev-Popov ghosts C, 
the antighosts C and the Lagrange multipliers B, and choose gauge-fixing conditions to make the 
functional integral perturbatively well defined. In addition, to implement the renormalization of 
divergences to all orders, study the gauge dependence and prove the Adler-Bardeen theorem, it is 
also convenient to introduce the sources K and use the Batalin-Vilkovisky formalism. The extra 
helds and the sources must be switched off at some point. In this section we explain how to dehne 
the physical quantities and show that they are gauge independent, under the sole assumption that 
the theory is AB nonanomalous, as in section 3. We work with convergent functionals, so we can 
set e = 0. We denote the e ^ 0 limits of T/j and the other functionals involved in our arguments 
by the same symbols used so far, since no confusion is expected to arise. 

First, we need to “un-gauge-hx” the theory, by switching off C, H and their sources Kq, Kb- 
This operation is regular inside the F functionals, once Feynman diagrams have been evaluated, 
but not inside the actions S and Sr, in the sense that if we un-gauge-hx the action, Feynman 
diagrams obviously become ill dehned. For this reason, some gauge dependence survives the un- 
gauge-hxing procedure. Besides un-gauge-hxing, we must switch off the sources K. The combined 
switch-off procedure allows us to dehne a physical F functional, identify its gauge symmetries, 
check that they close on-shell, and prove that no gauge dependence ahects the physical quantities. 

Since the gauge hxing is introduced by means of a canonical transformation, such as (2.5), 
when we vary the gauge-hxing parameters 6 = ^ we make a canonical transformation. Therefore, 
the equations (3.7) and (6.4) can be used to study the dependence of the physical quantities on 
the parameters 
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The information gathered so far is encoded in the key formulas 


(7.1) 

(7.2) 


(ri? 6 »,ri? 6 ») — 0 , 

j ■’ 

and is sufficient to achieve the goals of this section. We work on the e —>■ 0 limit of (6.4), rather 
than the one of (3.7), because everything we say starting from the former can be easily generalized 
to the other case. 


7.1 Quantum gauge algebra 


Formula (7.1) gives 


f Sr^Re Sl^Re 

J 5K^ 


f / SrSpie \ SiTug 

J \5K^ / 


l{iSRe,<^n) 


Si^Re 


(7.3) 


and tells us that F/jgi is invariant under the inhnitesimal (nonlocal) transformations 

Here and below w, w', etc., denote constant anticommuting parameters. Write <I>“ = {(/>*, (7“, (7“, il“} 
and Kct = to separate the classical helds (j/ and their sources from the 

extra fields and their sources. 

Observe that S is independent oi Kr and contains Kq only through the term — f This 

is also true after the canonical transformation (3.2), if we assume, for simplicity, that the functional 
(5(‘h, K') appearing in (3.2) is independent of Kq and Kr. Then Sg also satisfies {Sg, C) = B and 
{Sg,B) = 0. Moreover, the sources Kq and Kr cannot contribute to any nontrivial one-particle 
irreducible diagrams. Thus, after renormalization we still have {SRg,C) = B and {SRg,B) = 0, 
i.e. 5(7“ = wB°‘ and 5B^ = 0. 

Define 

tR{^) ^ rR,(d>,i^)le=B=i^=o> = S^'^\c=B=K=0- 

Observe that Tr((P) is independent of the ghosts C, because it has ghost number zero and after 
suppressing C and K no fields and/or sources of negative ghost numbers survive. For the same 
reason, which has ghost number equal to one, is linear in C. Clearly, 5C = SB = 0. Thus, 
when C, B and K are switched off, formula (7.3) turns into 


0 = 



j SlfRjcp) 

S(f>^ 


(7.4) 
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The terms proportional to SiThq/SC do not contribute to (7.4) because T is C independent. 
The terms proportional to SiTjig/5C and 5iTjig/6B disappear, because they multiply 6C and 5B, 
respectively. 

We call the “physical” T functional. The transformations encode the gauge sym¬ 

metry of Tij. Indeed, recall that 5(f>^ is linear in C and of course w. Replacing each ghost C 
with zu'A, where A(x) is a function having statistics opposite to the one of C, and dropping the 
products ww' after moving them to the left, we can define a symmetry transformation by 
the formula 






C-i^ro'A 


and prove, using equation (7.4), that T is invariant under this symmetry: 


SAtnicj)) = I = 0 . 


We call the quantum gauge transformations. To the lowest order in H they coincide with the 
starting gauge transformations, but at higher orders they are in general nonlocal functionals. We 
call the algebra of the transformations 6a quantum gauge algebra. 


7.2 Closure of the quantum gauge algebra 


Now we study the closure of the quantum gauge algebra. If we differentiate (7.1) with respect to 
K, we obtain 

(rK0,(5<h“) = o. 

Consider this equation in the case —)■ 6(p'^, then switch off C and B, and set AT = 0 at the 

end. Recalling that 6C = wB and 5B = 0, and observing that 6(p^ does not depend on and 
Kb, 'Are obtain 


f , f f 6r{d(j)^w') 5«f/j((()) 


(7.5) 


having multiplied to the left by w' and having defined 6'^°^ = 5<I>“ . The right-hand side 

of (7.5) is proportional to the (j>^ “T field equations”, which means that closure is achieved on 
shell. The left-hand side of (7.5) can be handled as follows. Since 64>^ and are linearly and 
quadratically proportional to the ghosts, respectively, we can write them in the form 


64r^ = w C^TM, 


SC^ = - 






where and T^_ are nonlocal functionals. Here the bar indices include the spacetime points where 
the corresponding fields are located and the summation over repeated bar indices understands the 
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integration over those spacetime points. Now, take formula (7.5) and replace (7“ with vj "+ 
A“ and being functions of the coordinates. The left-hand side of (7.5) is turned into 
ww'w"w"' times 

/ -1 

Finally, the whole formula (7.5) is equivalent 

[<5a,<5e]</>* = 5[aT]</>* + I u*^'(</>,A,S)^^^|^, (7.6) 

where 

[A,S]“ = j 

and A, S) are suitable functions. Formula (7.6) expresses the on shell closure of the quantum 
gauge algebra. 

The field transformations and the closure relations become clearer if we switch to a more 
explicit notation, where they read 

<5a(/>*(x) = J dSA^{y)T^[(t)]{x,y), [A,S]“(x) = J d^yd^zA'’{y)T.^{z)T^^[(l)]{x,y, z), 

T^[(p] and T^[</>] being (nonlocal) functionals that depend on two and three spacetime points, 
respectively. 


7.3 Gauge dependence of the physical T functional 

The last goal is to study the gauge dependence of Observe that the functional {Hro) 

that appears in formula (7.2) has ghost number equal to —1. Therefore, it must be proportional 
to the antighosts C and/or some sources K. This fact implies that the derivatives 
and 5i{Hji0)/SC^ are zero at C = K = 0. Moreover, {Hrq) does not depend on Kq and Kb, 
if the functional (5(d>, K') of (3.2) satisfies the same property, as we are assuming here. Setting 
C = B = K = d'm (7.2) we obtain 


de 


5fR(0) 

dx 


U*(0) 


dcj)^ 


(7.7) 


where 


u\4>) 


SriHm) 



C=B=K=0 


Formula (7.7) is the equation of gauge dependence satisfied by the physical functional f/j((/)). We 
can integrate it with the procedure described in subsection 6.2. The first term on the right-hand 
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side of (7.7) can be absorbed into redefinitions of the parameters A, while the second term can 
be absorbed into a change of field variables. We can do this for each gauge parameter 6, taking 
one at a time. We obtain that there exists redefinitions X{X',9) and a change of field variables 
(f>{(p',X',9) such that the transformed physical functional 

t'ji{<p',X')=fR{cPi<p',X’,9),XiX',9),9) 

is 9 independent. Setting 0 = 0 we get T'R{(f>',X') = A', 0), which in the end allows us to 

write 

tRi<f>,X,9)=tR{<p'i<p,X,9),X'{X,9),0). 

Since the entire gauge dependence is encoded into changes of field variables and redefinitions 
of parameters, it cannot affect the physical quantities contained in Tr((P). 

7.4 Unitarity 

In this subsection we prove (perturbative) unitarity, to emphasize why gauge independence is 
so crucial. For definiteness, we illustrate our arguments in Yang-Mills theories, but everything 
we say can be applied to quantum gravity, as well as any general gauge theory. We recall that 
perturbative unitarity is the statement that the identity SS^ = 1 holds diagrammatically, order 
by order in the perturbative expansion [32]. A necessary condition is that the free-held theory we 
perturb around propagates only physical degrees of freedom. A necessary and sufhcient condition 
is that when the identity SS^ = 1 is written as a cutting equation no unphysical degrees of 
freedom contribute to the cut propagators. 

There exists no gauge-hxing conditions where both unitarity and the locality of counterterms 
are manifest. If we want manifest unitarity, propagators must have only physical poles. This 
happens when we choose gauge-hxing functions of the Coulomb type, such as G{4>) = diAi, 
where i,j, ■ ■ ■ are space indices, inside the gauge fermion 'I'(d>) of (2.4). However, the locality of 
counterterms is not manifest in that gauge, since the Coulomb propagators contain denominators 
whose dominant terms (those that determine their ultraviolet behavior) do not depend on the 
energy (or do not depend on it in the correct way). Then, when we differentiate a Feynman 
diagram with respect to the energies of its external legs, the overall degree of divergence is not 
guaranteed to decrease, so we cannot prove the locality of counterterms in this way. Besides 
having a bad power-counting behavior at high energies, the propagators of the Coulomb gauge 
generate spurious divergences that are difficult to handle. 

To have a good power-counting behavior we need to equip the propagators with extra poles, 
some of which are unphysical. This is achieved for example by choosing the Lorenz gauge-fixing 
function G{4>) = in (2.4). The Fadeev-Popov ghosts then also have poles. The locality of 

counterterms is manifest, but unitarity is not. 
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The extra poles must cancel somehow, but their mutual compensation is not evident. The 
best way to prove this compensation is to use the gauge independence of the physical amplitudes, 
which allows us to switch back and forth between gauge-hxing conditions of the Lorentz type 
and gauge-fixing conditions of the Coulomb type. The former make the locality of counterterms 
manifest and hide unitarity, while the latter make unitarity manifest and hide the locality of 
counterterms. 

For example, choose the gauge fermion 

= y - diAt + , (7.8) 


which contains two gauge-hxing parameters, ^ and This functional interpolates between the 
Lorenz gauge ((" = 1) and the Coulomb gauge (C = 0). After integrating B out, the propagators 
of the gauge helds are 


{Ao{k)Ao{-k)), = -^{(E^ - -e), 

{Ai{k)Aj{-k))f^ = ^2 _ ^2 ^ 


{A,ik)Ao{-k)), = ^^{C-OEk 


k^P{k) 


(7.9) 


where k"^ = kiki and 


P{k) = C(C-E^ - ik^f - P(1 - 


while the ghost propagator is 

{C(k)C{-k)) = (7.10) 

We see that the propagators are well behaved, from the point of view of power counting, whenever 
7 ^ 0. They are not well behaved for C = 0; which is the Coulomb limit. The parameter C is a 
sort of cutoff that regulates the spurious divergences of the Coulomb gauge. Moreover, at (" = 0 
P{k) is equal to ^^(fc^)^ and only the physical poles survive. Instead, unphysical poles are present 
whenever C 7^ 0 - 

In the previous sections we have proved that the physical quantities are gauge independent. 
In particular, they are independent of ^ and C,. Thus, they are also unitary, and obey the locality 
of counterterms. We see that they are unitary by taking C = 0- We see that they obey the locality 
of counterterms by taking C 7^ 0- 

In the case of the standard model in flat space, we can easily generalize the proof of the 
Adler-Bardeen theorem given in ref. [6] to the family of gauge fermions (7.8), because they are all 
renormalizable. Then, the remarks of this subsection allow us to infer that the standard model 
in flat space is perturbatively unitary. 
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In ref. [7] a more general proof of the Adler-Bardeen theorem was given. It holds in a large class 
of nonrenormalizable theories, which includes the standard model coupled to quantum gravity. 
Combining the results of [7] with those of section 3, we can extend the validity of the Adler-Bardeen 
theorem to the most general local gauge fermions. In particular, using an analogue of (7.8), to 
switch between the Lorenz and Coulomb gauges of diffeomorphisms and Yang-Mills symmetries, 
we infer that the standard model coupled to quantum gravity is unitary as a perturbative quantum 
field theory. So are its extensions, as long as they satisfy the assumptions we have made. 

We stress again that gauge independence is crucial to reach these conclusions, since the Adler- 
Bardeen theorem per se ensures gauge invariance, but not gauge independence. 

8 Checks of high-order calculations based on gauge independence 

In this section we discuss how to use the results of this paper to check high-order calculations, 
under the assumptions of section 6. We have proved that 

Proposition 1 The beta functions of the physical parameters A may depend on the gauge param¬ 
eters f, but that dependence can always be reabsorbed into finite A redefinitions. 

This proposition also reminds us that there exists a class of subtraction schemes where the 
beta functions are gauge independent, in agreement with the general theorem proved in section 
3. If we are extremely lucky, the framework we choose to simplify high-order calculations might 
belong to that class. In ordinary situations, we may expect to be lucky only to the lowest orders, 
which may mean till three or four loops, or for special choices of the gauge fixing. However, we 
may not be able to identify the right framework in advance. Therefore, contrary to the usual lore, 
in general we cannot make checks of high-order calculations based on the assumption that A beta 
functions are completely gauge independent. 

Nevertheless, the beta functions cannot be gauge dependent in an arbitrary way, precisely 
because their gauge dependence must disappear in a suitable class of subtraction schemes. Thanks 
to this, a criterion to make checks of high-order calculations, based on gauge independence, still 
exists. It amounts to verify that every f dependence contained in the A beta functions can be 
cancelled by means of finite A redefinitions. In this section we show that the correct criterion, 
although less powerful than expected, is nontrivial and powerful enough. 

For definiteness, consider the standard model in flat space, and let A* collect the cou¬ 
pling, the squared gauge couplings, and the squared Yukawa couplings. The most general A beta 
functions have the form 

OO 

= ( 8 . 1 ) 

n=2 
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where XUi-'-in constants and the powers of h are inserted to emphasize the order of the loop 
expansion. The most general perturbative A redefinitions can be parametrized as 




( 8 . 2 ) 


n=2 


where are other constants. We have 


/?*' = X] + X] • • • Ai„_i ^ ^Afcl-fc^inAfci ' ' ' Afc^ 

n=2 n=2 m=2 

oo 

(8.3) 


n=2 


Proposition 1 ensures that the gauge dependence contained in the beta functions /3j can be ab¬ 
sorbed inside the redefinitions ( 8 . 2 ), that is to say there exist constants such that the 

couplings A( have gauge independent beta functions /3'^. Using this piece of information, we can 
determine which nontrivial checks of high-order calculations are available. 

The one-loop coefficients Xiihi cannot be changed, because they are scheme independent 
{Xi^i^i — Xiihi)- Therefore, they are also gauge independent. Comparing (8.1) and (8.3), we find 
that the other coefficients are related by the formula 


Xii---i„i Xil---ini + {n ^)^ij{ii---in-2Xin-liri}j j {h ■ ■-in-lXinjji + ' ' ' ) ( 8 - 4 ) 

where the dots stand for contributions involving with k < n. We can define an iterative 

procedure to determine by assuming that the constants with k < n are known, and 

requiring that Xi^...i^i be gauge independent. 

Now, if the number of couplings A is N, the tensors have = iV(X|^) independent 

components [33], while the tensors have components, where I is the number of 

loops. For N = 1 (that is to say a single coupling A) and £ > 2 it is always possible to absorb the 
gauge dependence into A redefinitions (as long as the one-loop coefficient x of tbe beta function 
does not vanish), because = Ci^£_i = 1. For £ = 2 it is not possible, because the second and 
third terms on the right-hand side of formula (8.4) cancel each other. Thus, two nontrivial checks 
are available for = 1, due to the gauge independence of the one-loop and two-loop coefficients 
of the beta function. 

For > 1 more nontrivial checks of high-order calculations based on gauge independence 
are available, because cjsf/ > Proposition 1 implies that the number of ^-independent 

components of the tensors Xii---ie+ii is obtained by modding out the redefinitions (8.2). Generically, 
this operation leaves 

cn/ — CN,e-i = N 
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£+1 



independent checks at i loops. This number is {i + N)/{N — 1) times less than the number we 
would obtain if the beta functions were completely gauge independent. Indeed, in that case we 
would have independent checks at I loops, which is equal to the number of constants 

So far, the beta functions of the standard model have been calculated to three loops [34] 
and the results are fully independent of the gauge-fixing parameters. Presumably, the convenient 
gauge-fixing functions and the clever treatments of the matrix 75 used in refs. [34] project onto 
the class of subtraction schemes where the beta functions are already gauge independent, at least 
to the lowest orders. However, we may expect that this coincidence will stop, sooner or later. 
When that happens, we must be aware of the facts pointed out in this section. Moreover, we stress 
that in the proofs of properties to all orders, such as the proof of the Adler-Bardeen theorem in 
nonrenormalizable theories [7], it is often more convenient to use subtraction schemes that are 
less practical from the calculational point of view, but more convenient from the theoretical side. 
There, it is also important to keep in mind that the beta functions do not need to be gauge 
independent. 

9 Conclusions 

In this paper we have derived generalized Ward identities for potentially anomalous theories, and 
used them to study the problem of gauge independence. The new equations contain an extra 
term that is responsible for a number of interesting effects. We have renormalized the equations 
of gauge dependence and integrated them. The result is that every gauge dependence can be 
absorbed into a canonical transformation acting on the renormalized T functional, provided that 
the finite local counterterms are appropriately fine-tuned. RG invariance is preserved and, as 
expected, the physical quantities are gauge independent. Nevertheless, the beta functions of the 
couplings may in general depend on the gauge choice. Gauge independence is useful to switch 
back and forth between gauge conditions that exhibit perturbative unitarity and gauge conditions 
that exhibit a correct power-counting behavior and the locality of counterterms. 

In several cases, the Adler-Bardeen theorem ensures that the gauge anomalies cancel to all 
orders, when they are trivial at one loop. However, it is not sufficient, 'per se, to ensure that the 
physical quantities are independent of the gauge fixing. In this paper we have proved that, in the 
end, gauge invariance does imply the gauge independence of the physical quantities. Precisely, we 
have shown that it is possible to renormalize the theory and fine-tune its finite local counterterms 
so that the cancellation of gauge anomalies ensured by the Adler-Bardeen theorem is preserved 
for arbitrary values of the gauge parameters. 

Said differently, assume that the gauge anomalies vanish for some specific choices of the gauge 
parameters. Varying or turning on a gauge parameter is equivalent to making a canonical trans¬ 
formation. After the canonical transformation, it is always possible to re-renormalize the theory 
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and re-fine-tune its finite local counterterms to enforce the cancellation of gauge anomalies again. 
Moreover, the gauge dependence of the renormalized F functional is encoded into a convergent 
canonical transformation. The theorem proved in section 3 is very general, to the extent that we 
did not need to make particular assumptions about the gauge algebra or the properties of the 
theory under renormalization. In particular, it holds for renormalizable and nonrenormalizable, 
chiral and nonchiral, theories and for arbitrary composite fields. Once we know that the cancel¬ 
lation of gauge anomalies holds in the framework we prefer, we know that it holds in every other 
framework. 

One application of the theorem is to power-counting renormalizable chiral gauge theories 
gauge-fixed by means of a nonrenormalizable gauge fixing. It allows us to show that the param¬ 
eters of negative dimensions introduced by the gauge hxing do not propagate into the physical 
quantities. In other words, the theory remains renormalizable, although in a nonmanifest form. 
A second application is a crucial step in the proof of the Adler-Bardeen theorem for nonrenor- 
malizable theories elaborated in ref. [7]. 

It is often possible to prove the cancellation of gauge anomalies in a family of gauges. In that 
case, if the assumptions listed in section 6 hold, we do not need a new fine-tuning to enforce the 
cancellation of gauge anomalies after the variation of a gauge parameter. Then, the gauge depen¬ 
dence of the theory is encoded into a convergent canonical transformation on the renormalized 
r functional, combined with a finite redefinition of the parameters. This fact makes it apparent 
that in general the beta functions of the couplings may depend on the gauge fixing. We expect 
that high-order calculations of the beta functions in the standard model will exhibit, sooner or 
later, dependences of the type mentioned here. 

The gauge dependences of the beta functions can be eliminated by redefining the couplings in 
ad hoc ways. Thanks to this fact, gauge independence can still be used to make nontrivial checks 
of the calculations. 


Appendices 

A Useful formulas 


In this appendix we collect a few identities that are used in the paper. First, we recall that 

(r,r) = ((5,5)), (A.i) 

where S is any action (renormalized or not), T denotes the F functional associated with S and 


(U = I [dd>]Xexp (f5(d>,A:) +z J d-V, 


(A.2) 
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is the average defined by S, X being a local functional. Formula (A.l) can be proved by making 
the change of field variables (2.10) in the functional integral (2.8), and recalling that in any dimen¬ 
sional regularization the local perturbative changes of field variables have Jacobian determinants 
identically equal to one. For details on the derivation, see the appendices of refs. [6, 8]. 

If ^ is any parameter, we also have the formulas 


ac“\ac/’ 

{r,{X)) = {{S,X))+^-{{S,S)X)r, 


(A.3) 

(A.4) 


where X is an arbitrary local functional and {XY)-^ denotes the set of one-particle irreducible 
diagrams that have one X insertion, one Y insertion, and arbitrary and K external legs, Y 
being another local functional. Formula (A.3) follows from the definition of F as the Legendre 
transform of W. Formula (A.4) can be proved by making the change of field variables (2.10) in 
the average (A.2), and expressing the final result in terms of <I> and K. For details on this method, 
see the appendix of ref. [8]^. 

A simpler method to derive formula (A.4) is to deform the action S into S + Xa, where 
O' is a constant, consider the deformed version of formula (A.l) and take the first order of its 
expansion in powers of a. By (A.3), F is deformed into F -|- {X)a + O(cj^). Instead, the average 
(y) of a local functional Y is deformed into (Y) + i{YX)YO' + 0((T^). Indeed, the factor e*"^ 
appearing in the integrands of Z{ J,K) and Z{ J, K) (Y) [check (2.8) and (A.2)] is deformed into 
e*S(i -\-iXa + 0(cj^)). Moreover, the deformed average, considered as a functional of and K, is 
still a collection of one-particle irreducible diagrams. Thus, the first correction to (Y) is precisely 
i{YX)rcr- Taking Y = (5,5), we obtain ((5,5)) ^ ((5,5)) -|-i((5, 5)X)rO'+ ©(cr^), wherefrom 
(A.4) follows. 

If we subtract the equations (A.3) and (A.4) we also get 


^ - (r, (y)) = {^- - ^(<5,5)A^^ , (A.5) 

which is the starting point to derive the equations of gauge dependence. 

Another useful identity tells us that [16, 8], if —)■ ^',K' is a canonical transformation 

with generating functional F{^,K'), and Y{^,K) is a functional behaving as a scalar, i.e. such 
that F'(4>', AT') = Y{^,K), then 


dY' 


dY 




(A.6) 


where F^(<h, iF) = F^(<I>, iF'(<I>, AT)) and A<^(<I>,Ar') = dF/dC,. The field and source variables that 
are kept constant in the C, derivative of a functional are the natural field and source variables of 
that functional (that is to say <!>' and K' for Y', and K for T, and K' for F). 

^Note that we have switched from the Euclidean notation used in [8] to the Minkowskian notation used here. 
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B Renormalization of local bifunctionals 


In this appendix we show how to renormalize a generic local bifunctional, and then specialize to 
evanescent local bifunctionals. Given a theory with action S, assume that a local bifunctional S' 
has the form AB, where A and B are local functionals. Couple A and B to external (constant) 
sources Ha and hs, by deforming the action S into S = S — iHaA — ihsB. Then, renormalize the 
extended action S. The renormalized version of S has the form 

Sr = Sr — BiaAr — iJirBr — ihshACR + 0{h\) + 0{h%), 

where Ar and Br are the renormalized functionals A and B, respectively, Ha, hs are the renor¬ 
malized sources, and Cr is a local functional. Consider the T functional Tr associated with Sr. 
Differentiating it from the left-hand side with respect to hs and then Ha, and later setting Ka = 
Kb = 0, we find that the renormalized S' is equal to 3^r = ArBr -|- Cr. 

It is a known fact (see for example [21], chapter 13, or [6], section 6) that an evanescent local 
functional E can be renormalized so that its renormalized version Er satisfies {Er) = 0(e). This 
property extends to evanescent local bifunctionals in a straightforward way. However, we have to 
pay attention to some details. 

By writing and everywhere inside E, we can express each vertex of E 

in a factorized form where 6k denotes the evanescent part, made of tensors , possibly e 

factors and other structures that stay outside of the diagrams, while Tfc is a nonevanescent local 
functional and collects all the momenta. We then have E = Instead of considering 

the average {E), consider first the diagrams (T^) that contain one insertion of T^. Iterating in 
n = 0,1,..., let denote the (n -|- l)-loop divergent part of {7nk), where 

n 

p=i 

are the functionals renormalized up to and including n loops. By the locality of counterterms, 
each is local. Then, the functional En = is renormalized up to and including n 

loops, and satisfies 

{En) = {7nk) 6k = 0(e) + 0(r+i), (B.l) 

k 

because each (T„fc) is convergent up to Finally, the functional Er = E^^ satisfies 

{Er) = 0(e). 

(p) 

In the procedure just outlined we have subtracted away all sorts of contributions order 

by order. More generally, we do not need to subtract those that, once multiplied by 6k, give 
evanescent results. Indeed, collecting those evanescent local parts inside a local functional AE, 
anything we have said so far for E can be repeated for AE. We reach the conclusion that 
{Er) = 0(e) even if we “forget” to subtract any evanescent local parts. 
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Once we have renormalized E so that (Eji) is evanescent to all orders, we can apply the same 
procedure to the bifunctional Y = EB, where B is an arbitrary local functional. The outcome is 
that we can find a 0(h) local functional Tr, such that the local bifunctional Yr = Ej^Bji + Tr is 
renormalized and the average (Yr)p is evanescent to all orders. 

More precisely, we can iterate the renormalization of Y as follows. Write Y = EB = SkUk, 
where 11^ = T^i?. Let Bn denote the functional B renormalized up to and including n loops. 
Inductively assume that the n-loop renormalized U/. have the form Unk = ‘^nkBn + Onk^ where 
Qnk are local functionals. Define Yn = Ylk^kUnk = EnBn + Fn, where = Ylk^kOnk- Clearly, 
(^n)r = 0(e) + 0(hP‘^^), because each ('lX„fc)r is convergent up to 0{hP'^^). By the locality of 
counterterms, the (n + l)-loop contributions to (Unk)r are made of a local divergent part 

^nfcdiv ’ plus a generically nonlocal convergent part. Consequently, the (n + l)-loop contributions 
to (Yn}r are the sum of a local divergent part, a local nonevanescent part, plus a generically 
nonlocal evanescent part. If we define Un+ik = '^n+ikEn+i YSn+iki where Bn+i is the functional 
B renormalized up to and including n + 1 loops, and Qn+ik = Onk (i^n+ifc)r 

is convergent up to 0{h"’~^‘^), and so (l^+i)r = 0(e) + 0{h'^~^‘^), where = Ylk^k'^n+ik = 

En+iBn+i + Fn+i, and F„+i = SkQn+ik- The conclusion also holds if we “forget” to subtract 
any evanescent local parts of En+i and/or Fn+i- The subtraction can be iterated in n so that in 
the end (Yr)p is evanescent to all orders in h, where Yr = Y^o- 


C Integrating equation (3.33) 


In this appendix we integrate the equations (3.33) and (3.32). First, we recall how to integrate 
the simpler equation 


dX 

de 


{X,V), 


(C.l) 


for the functional X{^,K,9), given the functional V{^,K,9). Expanding in powers of 9, write 


OO 

Vi^,K,9) = Y,0"Vn{^,K). 

n=0 

We want to show that there exists a canonical transformation <h,iL —)■ with generating 

functional 

/ OO 

<!>“< +^rT„(cI>,iL'), (C.2) 

n=l 

such that 

Y'($', K') = X($(^>', K', 9), K(^', K', 9), 9) 

is independent of 9. 
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We can derive conditions on the unknown functionals 3“^ by applying formula (A.6), which 
relates the functional V of (C.l) to the canonical transformation S'. A sufficient condition to have 
dX'jdO = 0 is y = ^ 0 , where Iq = dl/dO. In other words, 

0 = f; r [Vn{^, K)-{n + 1) Jn+I (<h, K')] , K^ = K + Y^ r 

n=0 n=l 

The first equation can be solved for T^+i by working recursively in n. It is sufficient to express each 
14(<h, AT) as a functional of $ and K\ by using the second equation, and then set the coefficient 
of 0"' to zero. This proves that the desired canonical transformation (C.2) does exist. Clearly, 
coincides with A(<I>', A'', 0). Therefore, expressing everything by means of fields and 
sources without primes, we get 

A(ci>, K, 9) = X(<I>'(cl>, AT, 9),K'{^, AT, 9),0). 

(C.3) 

Define a new functional G and a map 

G(<h, K, 9) = f d0G(a>, K, 9), K, 9) = T d9 (Z,-, , 

Jo Jo 

where Zg = Z{^, K, 9) and Vg = V (<i>, K, 9). Observe that 

Then, equation (C.3) turns into equation 

BY °° 

— = (y,c), fory = y- 

n=0 

Note that the terms are at least 0(0"'^^). Using the result found above, the canonical 

transformation <h, A' —<h',A'' given by formula (C.2) is such that the transformed functional 

F'(^>', K') = y(^>($', K', 0), Ar($', AT', 9), 9) 

is 9 independent. Finally, if G = 0{u'^) for some expansion parameter u (which is e or h, when 
we apply this theorem in subsection 3.2) and V is regular in u, then the canonical transformation 
<h, AT —)■ K' is also regular in u, which implies 

y($($', K', 9),K{<^', K', 9), 9) = K') + 0{u^). 


Now, assume that a functional Y[^, K,9) satisfies 

BY 

^ = (r.r) + G. 

where V (<h, K, 9) and G(‘h, AT, 9) are two other functionals. 
£0 : Z —)■ where Z is a functional, as 


47 



Setting 0 = 0, we get 

F($', K\ 0) = Y'{^\ K') + OK). 

Hence, expressing everything by means of fields and sources without primes, 

y($, K, 9) = y($'(^>, K, k, e),o) + oK). 

In other words, the functional Y{^,K,6) still evolves by means of a canonical transformation, 
but only up to 0{u'^). 

In most applications, the functionals V and G of equation (C.3) may intrinsically depend on 
Y. For example, this happens when Y is some renormalized action (or the F functional associated 
with it) and V, G are (the averages of) some renormalized local functionals, calculated with 
that action. We can disentangle this difficulty by expanding each functional in powers of h and 
proceeding inductively in this expansion. Writing 

OO OO OO 

Y = Y, ^ = E 

n=0 n=0 n=0 

we obtain the equations 

dY 

- (y„, vb) = (C.4) 

k=0 

which have the same form as (C.3). The contributions Vk and G^ to V and G with k ^ n do not 
depend on Y^- For k = 0 this is obvious. For A: > 0 it is sufficient to observe that the vertices 
Yn of order of the renormalized action Y can only contribute to the one-particle irreducible 
diagrams associated with V and G that have n -|- 1 or more loops. Indeed, at least one additional 
loop must be closed to connect a vertex Yn with the insertions provided by V or G. When Y is 
the F functional and V, G are averages of local functionals, we can argue similarly. 

Now, assume that we have solved the equations (C.4) for n < h, and consider the equations 
(C.4) for n = n. The unknown is Yn^ while 14 and G^ with k ^ n are independent of it. Thus, 
equations (C.4) can be solved with the method explained above. We conclude that the procedure 
we have given to solve the equations (C.3) is well defined. 

D Standard model coupled to quantum gravity 

In this appendix we report some reference formulas for the standard model coupled to quantum 
gravity. The classical fields cj) contain the vielbein e^, the Yang-Mills gauge fields and the 
matter fields, where the indices a,b,... refer to the Yang-Mills gauge group (within which we 
include the Abelian subgroup) and a,b,... refer to the Lorentz group. The classical action Sc{4>) 
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is equal to the sum S'cSM + ^Sc, where 


5'cSM = / \/W\ 


-^{R + 2A,)-^F^,F-P^ + Lr 


and ASc collects the invariants generated by renormalization as counterterms, multiplied by 
independent parameters. Here, R is the Ricci curvature, g is the determinant of the metric 
tensor, Fj^p is the Yang-Mills field strength, Lm is the matter Lagrangian coupled to gravity, Ac 
is the cosmological constant, and k? = SttG, where G is Newton’s constant. 

The functional Sk of formula (2.2) reads 


SK = j{CPdpAl + A}dpCP - dpC^ - gr'^^AlC^)Kf + j (^C%C^ + 

+ j\c%e^ + eppCP + C%iM + J CP{dpCnK? + j+ CPdpC°^)Kg 

+ I (c%4;l - + K^ + I(c%^l;L - 

+ I (Cf^idp^) + K^-J - I BpK^ - I B,j,Kg, 


where V’L are left-handed fermions, ip are scalars, while T“ and T" are the anti-Hermitian matrices 
associated with their representations. The triplets C^-C°'-B^, and C^-Cp-Bp collect 

the ghosts, the antighosts and the Lagrange multipliers of Yang-Mills symmetry, local Lorentz 
symmetry and diffeomorphisms, respectively. It is easy to check that {Sk,Sk) = 0 in arbitrary 
D dimensions. 

Finally, the gauge fermion of formula (2.4) reads 

'I'(^) = J VW\C^ + J eCpp (^eP-g~^^dpdJp + ^B^^ + 

- 1V\9\Cp + ^g'^^gp,d,gp-^ - ^gP^Bp^ , 

where ^-nd are gauge-fixing parameters. 


E Comparison with manifestly nonanomalous theories 

We have mentioned that an unexpected consequence of our results is that in AB nonanomalous 
theories the beta functions of the couplings can depend on the gauge-fixing parameters. It is 
interesting to better understand why this does not happen in manifestly nonanomalous theories. 

We actually begin with nongauge theories, that is to say theories that have no gauge sym¬ 
metries. There the action S{^,K) does not even depend on the sources K and the canonical 
transformations are just arbitrary changes of field variables. 
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Denote the classical action by S{(j)), the renormalized action by Sii{4>) the renormalized F 
functional by Fr((/>). We assume that Sr and F/j are defined by subtracting away the divergences 
just as they come, in the minimal subtraction scheme. 

Consider a local, perturbative change of field variables 

^P\<P,9) = <P^ + (D{e) (E.l) 

for the classical action S. Let Sg{(j),9) denote the transformed classical action, 

Sg{cP,9) = Sii;icP,9)), 


which obviously satisfies 


where 


dSe 

89 


= / A(/>^ 


SiSe 


= 




(E.2) 


S9 5^^ 

Denote the renormalized Sq by Sro and the F functional associated with it by F/j^. 

We want to show that the change of field variables (E.l) on S is mapped onto a renormalized 
change of field variables on Sr and a nonlocal, convergent change of field variables on F/j. This 
property is encoded into the equations of gauge dependence, which now read 


dSRe 

89 


= / Ar(I)‘ 


^iSrb 


dVRe 

89 


— / 


SiTro 


(E.3) 


where A^(/>* is the renormalized version of the composite field (E.2). Equations (E.3) are just 
particular cases of equation (C.l), and can be integrated with the method explained in appendix 
C. So doing, it is straightforward to prove that the 9 dependences of both Srq and F/jg are encoded 
into pure changes of field variables, with no redefinitions of parameters. 

We point out that the first equations of formula (E.3) are highly nonlinear in Sr 0 , because 
A/j(/>*, being a renormalized composite field, intrinsically depends on Srq. Nevertheless, with 
the inductive procedure explained in appendix C we can disentangle this dependence. Similarly, 
the equations satisfied by F/jg contain the average on the right-hand side, which is also 

determined by Sro. The procedure to integrate the equations of F^^i is basically the same as the 
one for Srq and is again given in appendix C. 

Formulas (E.3) can be proved by induction, using the minimal subtraction scheme. Let Sn = 
So + 0(/i)xpoles and A„i;A* = A(p^ + O(^)xpoles denote the action and the composite field (E.2) 
renormalized up to and including n loops. Assume that 

I A„0*^ = (E.4) 
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Clearly, this assumption is satisfied for n = 0. Moreover, in the minimal subtraction scheme 
is made of pure poles. 

Differentiating the F functional F^, associated with Sn, with respect to 0, we get 


dVn _ /dSn\ 

de \ de / 


d^X ( /^n4>\x) 


SlSn 


+ {d^n)n- 


(E.5) 


Now, 


6lSn 


(^Sn + i 



—Jjix) — i 


Si 




exp (iSn + i 



Using this formula inside (E.5) we can drop the last term by integrating by parts, because when 
the derivative Si/S4S{x) acts on A„(/>*(x) it gives zero in dimensional regularization. Finally, we 
obtain 


dVr, 

d9 


= - / J^ + {dlnU = / 


Sj^ 

S<j)'^ 


+ (3^n)r 


(E.6) 


Since Sn and A„())* are renormalized up to and including n loops, the (n + l)-loop divergent parts 
^ndTv^ and V* of F„ and (^Anfp^)^ are local. Moreover, the divergent part of {‘3ln)n 

coincides with the part of because din starts from 0{h^~^^) and it is just made of 

poles. Thus, taking the divergent parts of formula (E.6) we get 


dr 


(n+l) 

ndiv 


= A 


(n+i)^i SiSe 


+ 


r p(n+l) 

ndiv _|_ O'Or) 


de j ^ Scf)^ J ^ de 

Subtracting the divergences just as they come, we define 


- / A, 


,i SlSn 


Scf)^ 




(E.7) 


c _ c _ r(”'+i) 

*Jn+l — *Jn i ndiv ’ 


An+l</>' = Anf - A 


(n+l) 
ndiv ' 


Clearly, the F functional F^+i associated with S'n+i is renormalized up to and including n+l 
loops. Using (E.7), we find 

^ ^ 

Thus, the inductive assumption (E.4) is promoted to the next order. The equations (E.3) 
follow by taking n = oo in (E.4) and (E.6). 

We see that in theories with no gauge symmetries a change of field variables on the classical 
action does not generate redefinitions of parameters in the renormalized F functional: the pa¬ 
rameters e introduced by the field redefinition do not propagate into the beta functions of the 
couplings. Moreover, we do not need to re-fine-tune the finite local counterterms. 

Another approach to these issues was given in refs. [35, 36], where the changes of field 
variables were mapped from the classical action to the renormalized action and the (renormalized) 
generating functionals Z, W and F, as well as a more general type of F functional, called master 
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functional. That approach also shows that a change of field variables does not affect the beta 
functions of the couplings, in the theories that have no gauge symmetries. 

Similar properties hold in manifestly nonanomalous gauge theories, where the equations 

g™ = {Sji 0 , Qro), = (r^e, {Qne)) (E.8) 

hold and can be integrated [8]. Again, the conclusion is that a canonical transformation acting on 
the classical action is converted into a renormalized canonical transformation acting on the renor¬ 
malized action, and a nonlocal, convergent canonical transformation acting on the renormalized 
r functional, with no effect on the beta functions of the couplings. Equations (E.3) can also be 
obtained by switching off the sources K in formulas (E.8). 

What “goes wrong” in AB nonanomalous theories, is that “small things”, that is to say evanes¬ 
cent terms 0(e), are around all the time, and can generate unexpected finite corrections by sim¬ 
plifying some divergences. For this reason, they force us to re-fine-tune the subtraction scheme at 
every, even minor, modification of the framework in which we formulate the theory. Yet, we have 
shown in the paper that we can put their effects under control and preserve the correct physical 
properties. 
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